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Clock synchronization is the backbone of applications such as high-accuracy satellite navigation, 
geolocation, space-based interferometry, and cryptographic communication systems. The high ac- 
curacy of synchronization needed over satellite-to-ground and satellite-to-satellite distances requires 
the use of general relativistic concepts. The role of geometrical optics and antenna phase center 
approximations are discussed in high accuracy work. The clock synchronization problem is explored 
from a general relativistic point of view, with emphasis on the local measurement process and the use 
of the tetrad formalism as the correct model of relativistic measurements. The treatment makes use 
of J. L. Synge's world function of space-time as a basic coordinate independent geometric concept. 
A metric is used for space-time in the vicinity of the Earth, where coordinate time is proper time on 
the geoid. The problem of satellite clock syntonization is analyzed by numerically integrating the 
geodesic equations of motion for low-Earth orbit (LEO), geosynchronous orbit (GEO), and highly 
elliptical orbit (HEO) satellites. Proper time minus coordinate time is computed for satellites in 
these orbital regimes. The frequency shift as a function of time is computed for a signal observed 
on the Earth's geoid from a LEO, GEO, and HEO satellite. Finally, the problem of geolocation in 
curved space-time is briefly explored using the world function formalism. 
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I. INTRODUCTION 

In the past several decades, there has been a dramatic improvement in two technology areas: atomic clocks 
and lasers for free-space optical communications. Future spacecraft are envisioned as communicating by free- 
space laser links. The synergy of technological development in atomic clocks and free-space lasers will lead 
to unprecedented advancements in high-accuracy space-time navigation P, digital communications, geoloca- 
tion 0, 0, IE IE EL surveillance using space-based interferometers ]E |^ , coherent distributed-aperture sensing at 
high frequencies ITU IT^ Il3 . ITEl \wi , and cryptographic communication systems. 

Accurate clock synchronization is the backbone of these systems. Consider the dependence of geolocation accuracy 
on clock synchronization. Consider geosynchronous satellites that receive a signal from an emitter of electromagnetic 
radiation located on the surface of the Earth. Assuming that the signals travel on the linc-of-sight, the geometry 
leads to a maximum signal time difference of arrival between two satellites that is 19.64 ms, see Appendix A. The 
information on the difference of ranges, Al, between the emitter and each satellite, is contained in the maximum time 
delay that is equal to 19.64 ms. If the clocks in two satellites are synchronized only to an accuracy of, say 10 ns, the 
order of magnitude in the position error of the emitter. Ax, is given by the fraction of the range difference: 

10ns . , , 

Ax a; - 3000m = 1.5nm (1) 

19.64ms ^ ' 

An improvement in the clock synchronization translates directly into an improvement in position accuracy. For 
example, an improvement in clock synchronization by three orders of magnitude will produce a position accuracy on 
the order of 3 m. 

Applications such as a space-based interferometer, can have even more stringent requirements on clock synchro- 
nization. For example, a multi-satellite space-based interferometer (distributed aperture system) that operates at a 
wavelength A will require accurate determination of relative satellite positions (nodes of the interferometer) to better 
than Aa; = A. For practical purposes, say we will need Ax = A/10, which translates to a clock synchronization 
requirement (A/lOc), where c is the speed of light. As an example of the stringent requirements on time synchroniza- 
tion, consider operating at 22 GHz, which is a frequency of interest to radio astronomers |l^. At this frequency, the 
position of the satellite nodes must be known accurately to A/ (lOc) =1.4 mm and time synchronization to 4.5 ps. See 
Table m for a range of values corresponding to different frequencies. For operating at optical wavelengths of 500 nm, 
the required position must be known to 50 nm and time must be synchronized to 0.16 fs. These numbers challenge 
and exceed the realm of possible time synchronization accuracy that is available today. However, the synergy between 
accurate clocks and optical free-space communication is expected to continue, so that hardware may soon support 
such stringent time synchronization requirements. New time synchronization schemes will then be required. 

The emerging field of quantum information and quantum computation |T9l l20l | has potential to produce new ultra- 
precise clock synchronization protocols. In fact, several clock synchronization schemes have been proposed based 
on quantum mechanical concepts |^ I2E HE IS HE HE HE HE HE 113- However, most of these schemes have 
neglected real features of the clock synchronization problem that are essential for real-life applications: the clocks to 
be synchronized are in relative motion, and at varying gravitational potentials. The fact that clocks are affected by 
their motion and by the gravitational potential are basic concepts that have their origin in Einstein's special and 
general relativity theory. The examples of the required accuracy of clock synchronization cited above, together with 
the magnitudes of relativistic effects on satellites (see Section IV, subsection E), show that the required accuracies can 
only be met by theories that take into account the effect of gravitational potential differences and relative motion. If 
synchronization of clocks is to be achieved using quantum information concepts, then certain features, such as relative 
motion of clocks and effects of gravitational potential, must be incorporated in the quantum information approach to 
clock synchronization. 

This article addresses the theoretical problem of clock synchronizing and syntonization (making two clocks run at 
the same rate), or correlating time on clocks on satellite platforms that are orbiting Earth, or that are near-Earth. 
However, Lorentz transformations between two systems of coordinates in relative motion show that space and time 
are really interwoven. The problem of clock synchronization is really part of the more general problem of navigation 
in space-time Ij 2]. For example, in the Global Positioning System (GPS), a user's receiver determines three spatial 
coordinates as well as time 2j. Therefore, in this article we will focus on the complete problem of navigation in 
space-time from the point of view of curved space-time, such as is invoked in general relativity. The features that 
relativity deals with, motion of clocks and effect of gravitational potential, are also features that must also be included 
in any classical or quantum theory of space-time navigation. We use the fact that space-time is described by a four 
dimensional metric, but for the most part we do not explicitly use the field equations of general relativity. In this 
sense our discussion is not restricted to general relativity. 

In this article, several considerations are stressed in the space-time navigation problem. First, space-time navigation 
is based on real measurements made by real physical devices. Real measurements are (spatially and temporally) local 
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TABLE I: Frequency /, wavelength A, and A/(10c) time synchronization is shown for various operating regimes of an interfer- 
ometer. 



/ A A/ (10c) time synchronization 

1500 kHz 200 m 66ns 

20 MHz 15 m 5 ns 

120 MHz 2.5 m 0.83 ns 

20 GHz 1.5 cm 5.0 ps 

60 GHz 5 mm 1.6 ps 

30 THz 10 m 3.3 fs 

500 THz 500 mn 0.16 fs 



quantities that are invariants under change of space-time coordinates (see Section VI) . Historically, due to a lack of 
accurate measurements over large distances (such as spacecraft to ground) measurement theory has not played a large 
role in relativity theory. On the other hand, measurements are at the core of quantum theory. Perhaps this point 
of intersection between relativity and quantum theory will help clarify how to augment quantum information theory 
with relativistic ideas. 



II. HARDWARE TIME, PROPER TIME, AND COORDINATE TIME 

A clock is a physical device consisting of an oscillator running at some angular frequency w, and a counter that 
counts the cycles. The period of the oscillator, T = 27r/a;, is calibrated to some standard oscillator. The counter 
simply counts the cycles of the oscillator. Since some epoch, or the event at which the count started, we say that a 
quantity of time equal to NT has elapsed, if N cycles have been counted. 

In this article I distinguish between three types of time: hardware time r*, proper time t, and coordinate time t. 
Hardware time is associated with a real physical device that keeps time, which I call a hardware clock. Specifically, 
hardware time t* — NT is the time kept by a hardware clock, and is given in terms of the number of cycles N counted 
by the device. Two hardware clocks will differ in the elapsed time that they indicate between two events, because 
no two devices are exactly the same. Furthermore, heating, cooling, and vibration typically affects real devices, and 
consequently the value of the elapsed hardware time registered on a real hardware clock can vary for these reasons. 

Proper time is an idealized time interval occurring in the theory of relativity. We imagine that there exists an 
ideal clock (oscillator plus counter) that is unaffected by temperature or vibration. However, based on Einstein's 
theory of relativity ^J, the ideal clock is affected by gravitational fields, acceleration and velocities. According to 
Einstein's general theory of relativity, gravitational fields, acceleration and velocities affect all physical processes, and 
hence these effects are associated with the geometry of space and time. A basic tenet of the theory is that between 
any two events that are infinitesimally separated in space-time by dx^, i — 0,1, 2, 3, there exists an invariant quantity 
ds called the space-time interval js^l 

ds'^ = —gijdx^dx-' (2) 

where gij is the metric of the 4-dimensional space-time. In general relativity, if the two events are time-like, then, 
ds^ > 0, and the events are "separated farther in time than in space" . We interpret ds — cdr, where dr is the proper 
time elapsed on an ideal clock that moves between these two events. 

The definition of an ideal clock is one that keeps proper time intervals. More specifically, if a clock moves from 
point Pi to point P2 on a 4-dimensional world line given by coordinates a;*(u), i = 0, 1, 2, 3, for ui < u < U2, where u 
is some parameter such that Pi — {x^{ui)} and P2 — {a;*(u2)}, then the proper time interval At between these two 
events is given by the path integral in the space-time: 

I dx' dxi , j dx' dxi , 

where the second integral has been parametrized by the coordinate time length = ct, where t has units of time. 
From Eq. ©, it appears that the proper time interval At depends on the metric of the space-time, gij, and on the 
4- velocity of the clock, dx^ /ds. But in fact, At is a geometric quantity that depends on the path in space-time, and 
is independent of the coordinates used to compute At. For the purpose of applications, Eq. (PJ provides a functional 
relation between the proper time interval measured on an ideal clock. At, and the path which the clock traversed in 
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4-dimensional space-time. Under certain conditions. Eq. (j^J can provide a relation between elapsed proper time At 
and coordinate time interval, Ax^ = X2 — x^, 

Ar^hixlxl) (4) 

which depends on the coordinate time x\ and x§ of events Pi and P21 respectively. However, in general Ar depends 
on the whole path x^{u), and not just on the end points of the path. 

A good, real clock, which we will call a hardware clock, is believed to provide an approximate measure of elapsed 
proper time. Therefore, for a given real hardware clock, there is some functional relation between the elapsed proper 
time, Ar, and hardware time, Ar*: 

Ar* = /2(Ar) (5) 

where the function /2() has some stochastic aspects and is hardware dependent. 

A peculiarity of the theory of relativity is that the coordinate time, as well as the other coordinates x", 
a = 1,2,3, are not directly measureable quantities [sSl- The coordinates are only mathematical constructs, and 
cannot be measured directly. Events occurring in space-time have coordinates associated with them. Two events 
are labelled with coordinates Pi — {x\} and P2 = {x\] even though these quantities are not directly measureable. 
However, the theory allows us to relate the difference of coordinate time, x'^—x\^ between two events, to the proper time 
interval, by the path integral in Eq. Of course, as mentioned previously, neither of these quantities are measureable 
directly, instead, only hardware time intervals, Ar*, are measured directly from hardware clocks. Everything else is 
calculated. 

The relations in Eq. ©-(O permit, under some circumstances, the relation of measured hardware times Ar* to 
coordinates times Ax", which enter into the theory or relativity. 



A. Synchronization versus Syntonization 

The relation between proper time Ar and coordinate time is such that they may "run at different rates" . This 
is clearly the case when the metric of the space-time is not a constant over the integration path. For example, 
a good clock may measure proper time intervals accurately, however, because of Eq. this clock will run at a 
rate that differs from coordinate time because dr/dx^ ^ constant. While coordinate time is a global coordinate 
quantity, valid (almost) everywhere in the space-time, the proper time interval depends on the world line of the clock. 
For example, for an ideal clock that is stationary (has constant spatial coordinates x", a = 1,2,3) with world line 

= (x", x^, x^, x'^), for Xq < x° < x" +(ix°, the proper time interval is dr = ^J—gaodx^. Therefore, the rate of proper 
time, dr/dx'^, depends on position thorough 500- 

Consider now two ideal clocks at the same location and assume that these two clocks are synchronized to read the 
same starting time at some epoch, or starting event. Next move the clocks apart (hence they travel on different world 
lines) and then bring them together once again to a common location. On comparing the times on these clocks, we 
find that different amounts of proper time have elapsed on each of them. In other words, proper time ran at different 
rates for each of the clocks. We say that the two clocks were not syntonized (i.e., they did not run at the same rate) 
since, when they were brought back together a different amount of proper time elapsed on each clock. The rate of 
each clock can be compared at any instant to the underlying coordinate time (which is a globally defined quantity), 
by using Eq. (|3Jl, and in this way the time on one clock can be compared to the time on the other clock. In Section 
IX, we will find that the most serious problem for practical satellite applications is the syntonization of clocks. 



III. CHOICE OF A PHYSICAL THEORY 



A theory must be chosen as the basis for navigation and clock synchronization. The theory must deal with two 
areas: the space-time in which all events occur, and the propagation of electromagnetic fields in this space-time. 
While these two areas are coupled in the sense that electromagnetic fields are sources for the curvature of space-time 
(within general relativity), we will take these two areas separately. More specifically, we will neglect the (very small) 
effect of the electromagnetic field creating space-time curvature. We will assume that all curvature is caused by the 
presence of mass. For applications in the vicinity of the Earth, this is an excellent approximation. 

The most common theory to choose is general relativity, as developed by Einstein [sil [s^ • As has already been 
discussed, this theory includes the effects of motion on clocks and also the effects of gravitational potential on clocks. 
In fact, there are several variants of relativistic theories of gravitation. However, Einstein's general theory of relativity 
has so far passed all physical tests |35ll36j |. and we shall subscribe to it in this report. Einstein's general theory of 
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relativity can be divided into two parts: first, the theory assumes that there exists a space-time metric gij, given by 
Eq. Q, which relates proper time dr — ds/cXo the space-time geometric interval ds, and coordinate difference dx^ 
between two events. This idea is quite general, and is really an embodiment of the principle of equivalence [Sfills^. 
The principle of equivalence essentially states that gravitational mass (in Newton's universal law of gravitation) and 
inertial mass (in Newton's second law of motion) are the same quantity. Roughly speaking, the equivalence principle 
says that two small test bodies of different mass will fall along the same geodesic path, i.e., the same distance in the 
same time. This principle has been extensively tested [S^H^, and is the basis of essentially all metric-based theories 
of gravity, because they are based on geometric ideas embodied in Eq. Q . Physicists today generally agree that a 
theory of gravity should be a metric theory in a pseudo-Riemannian curved space-time with metric of the form given 
byEq. 0. 

The second part of Einstein's general relativity theory consists of the field equations |3lL 13^ 



G^^ = -kT*^ (6) 

where k = SttG/c^, where G is Newton's gravitational constant, and c is the speed of light in vacuum. The field 
Eqs. © relate the matter distribution, given by the stress energy tensor T'^ , to the effect that this matter has on 
space-time via the Einstein tensor. 



Gij - Rij - T^gijR (7) 



where the Ricci tensor, Rij = i?^^^,, is related to the Riemann tensor 

T^i j^i pi 1 pa pi pa pi fo\ 

jkm jm,k jk^ni ~^ ji^^ afc jk am \ / 



and the affine connection 



1 ^l 



r}fe = i9ji,k + gki.j - gjkj) (9) 

is related to the metric gij. So the field Eqs. l|njl are a set of equations for the components of the metric tensor field 
gij. In Eq. ordinary partial derivatives with respect to the coordinates are indicated by commas. 

The field equations of Einstein, given in Eq. ((HI, have only been solved and tested in a limited number of cases. 
Consequently, the field equations are on a less-firm footing than the equivalence principle. Fortunately, most of the 
applications of navigation and clock synchronization in a curved space-time rely only on the fact that space-time is 
a metric theory, given by Eq. |(5J). Therefore, our conclusions below regarding navigation and clock synchronization 
transcend general relativity. Specifically, our conclusions are based on the assumed-correctness of the equivalence 
principle. 



A. Electromagnetic Waves 

The electromagnetic field plays a central role in experiments and applications. In technology applications, all 
information is currently carried by travelling electromagnetic fields. Inter-satellite links, and ground to satellite links 
are all done using electromagnetic radiation fields. Time transfer between stations on the ground and satellites is 
done by electromagnetic fields. Since electromagnetic fields paly such a key role, in this subsection we briefly outline 
the equations governing the propagation of electromagnetic waves, namely the Maxwell equations in flat space-time, 
and their generalization to curved space-time. 

Electromagnetic fields in a medium, such as air, a dielectric, or a magnet, are described by the macroscopic Maxwell's 
equations, and in SI miits, in flat space-time, these equations take the form 

divD = p (10) 

divB = (11) 
dT) 

curlH J + — (12) 

curlE = (13) 
at 

where E and B are the electric field and magnetic induction (or magnetic field), respectively, and D and H are the 
electric displacement field and magnetic intensity, respectively. In a medium, the fields E and D, and the fields B 
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and H are related by constitutive relations. In a vacuum, these fields are simply related by: 

D = eoE (14) 
B = MoH (15) 

where Cq and /io are the permitivitty and permeability of vacuum, respectively. In an isotropic (but not necessarily 
homogeneous) medium such as the Earth's atmosphere, the constitutive equations are 

D = eE (16) 
B ^ fiH (17) 

where e and /x are the permitivitty and permeability of the medium, respectively. 

In a curved space-time, the physics of electromagnetic wave propagation has been less explored |3^l8^l3^l8^l4oLl4ll |. 

However, it is known that the gravitational field scatters and diffracts electromagnetic waves, and that the plane of 
polarization of an electromagnetic wave is rotated as the wave propagates through a gravitational field. In general, 
a gravitational field affects electromagnetic wave propagation similarly to a dispersive medium For weak fields, 
such as exist in the vicinity of the earth, these gravitational effects are smaller than the dispersive effects due to 
the atmosphere (at low altitude). The general equations that govern electromagnetic wave phenomena in curved 
space-time, in the presence of a dielectric are given by 40j 

F^J,k + FJk,^ + Fk^,J = (18) 

and 

H% = r (19) 

where Fij and H"^^ are two antisymmetric tensor fields, the comma indicates partial differentiation with respect to the 
coordinates and the semicolon indicates covariant differentiation with respect to the coordinates. All Roman indices 
take values i = 0, 1, 2, 3. The two fields are related by a constitutive relation 

ik 2 I ^ \ ^ia „,kb 



^^W^^c'\--j f^'^^'Fab (20) 

The constitutive relation in Eq (|20|l assumes that the medium is isotropic, but not necessarily homogeneous, so the 
permittivity e and permeability /i are both functions of position. In Eq H20(l . we have used the definitions of the 
contravariant components of the metric tensor g^^ and the effective metric 

-f'^ ^ g'^ - {n^ ~l)u'u^ (21) 

where n is the scalar index of refraction given by 

= ^ (22) 

where g = det and 7 = det7y, is the local 4- velocity of the medium in our system of coordinates, and J* is 
the 4-current density. In the proper frame of reference of the medium (where the medium is at rest), the field tensors 
take the simple forms: 

and 

= _F7 n w ' (24) 
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the current density is J' ~ {poC,J°'), where po is the proper charge density and J" is the current. In this proper 
frame of reference of the medium, Eqs. H18(l - (|19|l reduce to Eqs ()1() |I - H13|I . and the consitutive Eq. H2()|l reduces to 
relations in Eq. lfTC|l - lfT7|) . 
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If the medium is vacuum, then n = 1 and 7*-' = g*-', and the tensors Fik and Hik are not independent, instead they 
differ only by trivial constants of the vacuum. 

Equations lP)l - (P|l govern the propagation of electromagnetic fields in the presence of a medium in a gravitational 
field. In general, in the vicinity of the Earth, the dispersive effects of the medium are larger than those of the 
gravitational field. Due to the complexity of these e quat ions, they have not been explored in detail. Only several 
treatments have been attempted, see for example Refs. [HESEl. The Eqs. (CHIl-CSl) or |ini)~(IT3I), form the basis for 
applications such as time transfer, clock synchronization and communication. These equations are valid in a system 
of coordinates (frame of reference) that is in arbitrary motion. In particular, these equations describe propagation of 
electromagnetic waves, and the reception and transmission properties of antennas in the radio portion of the spectrum. 

B. The Geometrical Optics Approximation 

When discussing precise measurements, it is important to state precisely the theory and approximations used. Up 
to this point in time, the geometric optics approximation has been universally used in discussions of time transfer, 
usually without stating its use, and without examining the limitation of this approximation. Below, we state the 
geometric optics approximation and how it enters into time transfer ideas. 

The geometric optics approximation consists of the assumption that the wavelength of the travelling electromagnetic 
wave. A, is much smaller than the linear dimension I of all objects (length scales) in the physical problem under 
consideration [s^: 

X«l (25) 

In physical terms, the limit of short wavelength waves (geometrical optics limit) corresponds to waves that travel 
along straight lines, so that diffraction (e.g., bending of waves around boundary edges) is absent. As an example of 
geometric optics at work, consider the visible shadows cast on the ground by objects in the path of the light from the 
sun to ground. The light travels at approximately straight lines. However, if one looks very closely near the edge of 
the shadow, the boundary between light and dark areas is not sharp, and this is where geometrical optics shows its 
limitation-there is diffraction, or bending of the rays around sharp edges of an opaque material. 

In the literature, statements are often made in the context of special relativity theory that in flat space-time "light 
travels along a straight line" , or in curved space-time, that "light travels along a geodesic" . Both of these statements 
are true only within the context of the geometrical optics limit ji^, 14^1 143 . l45l | . In recent years, there has been limited 
work to explore the nature of travell ing electromagnetic waves in a curved space-time, going beyond the geometrical 
optics approximation [ssL |3^, "4^ 0, 0, |4^ . The gravitational field creates complex effects such as diffraction of 
the electromagnetic wave and rotating its polarization. 

C. Signal Detection and Use of Antenna Phase Center 

In precision measurements, it is important to have a clear concept of the point from which radiation emanates and 
the point at which the radiation is detected. Usually, such a discussion makes implicit use of the geometric optics 
approximation. 

One place where the geometric optics approximation is not well-satisfied is for real (radio frequency) antennas, 
because antennas are efficient at receiving and transmitting radiation at wavelengths that are comparable to the 
antenna size, so Eq. H25|) is not well-satisfied. The desire to continue to use the (very convenient) geometric optics 
approximation forces us to invoke the concept of antenna phase center in precise time transfer or navigation applica- 
tions. We then imagine that there is a unique emission point, from which radiation emanates, and a unique reception 
point, where the radiation is absorbed. 

The antenna phase center is defined as the apparent point from which radiation emanates (or is absorbed). In the 
far-field radiation region, for one vector component of the electric field of antenna radiation, and for one polarization 
and one frequency to, the electric field can be expressed as 

E = ui;(6',0)e'^(''''^)— e"*"* (26) 
r 

where the vector u is a real polarization unit vector, E{9, 0) is the (real) electric field amplitude, ip{d, (j)) is the (real) 
phase, r is the distance from the antenna, and i — \/— T. If a point can be found such that ip{9,(j)) is independent 
of 6 and (f>, i.e., independent of the direction of propagation, then this point is the antenna phase center. For most 
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practical antennas, no such point exists |46l l47j . The reception of an antenna is related to its transmission properties, 
with the same phase center, by reciprocity relations [i^lH]. 



where Re takes the real part of a complex expression, E is the electric radiation field at the receiving antenna, and 
h(n) is the receiving antenna vector effective length (sometimes called the effective height), which is a complex vector 
quantity that characterizes the electromagnetic wave phase relationship of the antenna in receive and transmit mode. 
The receive and transmit modes are related by reciprocity relations [43, HI, HE Hyl 

Despite the limitations of the concept of antenna "phase center" , this concept is routinely invoked in practice in 
real antenna systems. The practical complication is that the antenna "phase center" is not a fixed point, but its 
position depends on electromagnetic wave frequency, w, polarization vector u, and direction of travel with respect to 
the antenna (both in receive and transmit modes). In other words the phase center position is not fixed, instead it 
varies with lo and u over some range of coordinate values Aa; that is on the order of a wavelength of the radiation, 
Aa; ~ A. In other words, the point at which an antenna receives a signal is only precisely defined (in geometrical 
optics) to within a distance Aa; ~ A. Consequently, when one antenna transmits radiation and another antenna 
receives radiation, the effective distance between these two antennas can vary with frequency and polarization of the 
radiation, and with the relative orientation of the antennas. When the relative orientation of one satellite changes with 
respect to another satellite due to their relative motion, their effective separation changes due to change of relative 
orientation of their antennas, in addition to a real change of distance between them. Clearly, the concept of a single 
emission point and a single reception point is fuzzy on the order of the scale of a wavelength at both transmit and 
receive end points. 

In conclusion, the apparent distance between antennas changes with frequency, polarization, and orientation. In a 
case where precise navigation is to be carried out using radio frequencies, these effects must be taken into account, so 
that an accuracy of better than one wavelength of the radiation can be achieved. 

As an alternative to using radio frequencies, we can transition to satellites using optical frequencies, which have 
considerably shorter wavelengths. The transition to optical free-space communications in satellites, with smaller 
wavelengths on the order of 10 - 60 nm, has the advantage that we can invoke the geometric optics approximation, 
and suffer Aa; ~ A errors that are much smaller, because of the much smaller optical wavelengths. For example, for 
optical wavelengths in the range 10 nm to 60 nm, position errors are comparable to the wavelength, and this should 
be compared with radio frequencies, of say, 1 MHz to 10 GHz, with wavelength range of 3 cm to 300 m. 

Actually, at optical frequencies, Eq. H27() is not the physical mechanism that is responsible for detecting electro- 
magnetic fields. Instead, detectors of electromagnetic radiation work either as bolometers or quantum detectors. 
Bolometric detectors are based on the pyroelectric effect, which produces a change of dielectric polarization with 
increase of temperature, due to absorption of electromagnetic radiation |5lj |. The polarization change is detected 
electrically. 

The second common mode of detecting optical radiation is based on a quantum mechanical effect. For example, in 
a semiconductor material a photon (quantum) of the electromagnetic field is absorbed by an electron, and the electron 
makes a transition from a valence band quantum state to a conduction band state. The electron in the conduction 
band is then detected electrically. In either case, the minimum volume that is needed for detection is roughly of the 
dimensions of a wavelength of the optical radiation, so similar criteria apply (to that of radio frequency) for accuracy 
of the point of absorption. 

As discussed in the Introduction, we believe that future satellite systems will have optical links. Such satellites may 
have an interferometer that operates in the radio frequency portion of the spectrum, but the navigation (positioning) 
and time synchronization will be done optically. Consequently, the geometric approximation will be useful because 
of the (relatively) small optical wavelengths compared to radio frequency wavelengths. Therefore, we will freely use 
theory that assumes that electromagnetic waves travel on geodesies in curved space-time, which is the geometric optics 
approximation. 



The geometric optics approximation is valid for the applications addressed in this report. Within the geometric 
optics approximation we can say that electromagnetic waves travel on geodesies in space-time. This is a purely 
geometric statement: the emission and reception events of a signal are connected by a null geodesic. Geodesies in the 
space are determined by the metric of the space-time. 




V = Re{ h(n) • E } 



(27) 



IV. WORLD FUNCTION OF SPACE-TIME 
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A useful quantity to deal with measurements in space-time is the world function Q. The world function is simply one- 
half the square of the space-time interval (see Eq. Q), measured along the geodesic connecting two points [L..2..3L. 32| . 

The world function was initially introduced into tensor calculus by Ruse , Synge 0] , Yano and Muto [^3 j 

and Schouten |56j |. It was further developed and extensively used by Synge in applications to problems dealing 
with measurement theory in general relativity [sij . The world function has generally received little attention in the 
literature, so we provide a detailed definition here. Consider two points, Pi and P2, in a general space-time, connected 
by a unique geodesic path T given by x*(m), where ui < u < U2. A geodesic is defined by a class of special parameters 
u' that are related to one another by linear transformations u' — au + 6, where a and b are constants. Here, m is a 
particular parameter from the class of special parameters that define the geodesic F, and x'^{u) satisfy the geodesic 
equations 



■'^ du du 



The world function between Pi and P2 is defined as the integral along F: 



niPl,P2) 



2(^2 



Ml , 



dx' dx^ 



(28) 



(29) 



The value of the world function has a geometric meaning: it is one-half the square of the space-time distance between 
points Pi and P2. Its value depends only on the eight coordinates of the points Pi and P2- The value of the world 
function in Eq. H29II is independent of the particular special parameter u in the sense that under a transformation 
from one special parameter u to another, u' , given by m = au' + b, with x^{u) = x^(u(u')), the world function definition 
in Eq. (|29|l has the same form (with u replaced by u'). 

The world function is a two-point invariant in the sense that it is invariant under independent transformation of 
coordinates at Pi and at P2- Consequently, the world function characterizes the geometry of the space-time. For a 
given space-time, the world function between points Pi and P2 has the same value independent of the coordinates 
that are used. As a simple example of the world function for Minkowski space-time, consider 



1 



fl{x\,X2) — - r]ij Ax^ Ax^ 



(30) 



where rjij is the Minkowski metric with only non-zero diagonal components (— 1,+1,+1,+1), and Ax^ — {x\ — x\)^ 
1 = 0,1, 2, 3, where x\ and x\ are the coordinates of points Pi and respectively. Up to a sign, the world function 
gives one-half the square of the geometric measure (the interval) in space-time. Calculations of the world function for 
specific space-times can be found in Refs. 0, 0, |^ 113, IH, |5^ and apphcation to Fermi coordinates in Svnge|3l| and 
Gambi et al. "go]. 

The world function for the Schwarzschild metric, linearized in small parameter GM/c^r^ is given by 0: 



^^'qijAx'Ax^ 



GM 



|X2 - Xi| 



Xi I (cos Oi — cos 6*2) 



X2 - Xi 



log 



tan(^) 
tan(f) 



(31) 



where cAt 



C]', and Qi and Q2 are defined by 



cos d„ 



Xa ■ (X2 - Xi) 
|Xa||x2 - Xi| 



a = 1,2 



(32) 



The world function in Eq. (|31|) assumes a spherical Earth (J2 — 0), see Eq. (|69|) . For most applications dealing with 
delays of electromagnetic wave propagation, this approximation is sufficient. A more accurate calculation of the world 
function for Schwarzschild metric is given by Buchdahl and Warner 59] . Including the small effect of the oblateness 
of the Earth will require additional computations and is left for future work. For a pedagogic discussion of applying 
the world function to space-time navigation, see Ref. . 



A. Navigation in Curved Space-time 



A satellite must be localized with respect to some system of coordinates. The satellite can receive signals from four 
electromagnetic beacons and use the measurements to solve for its position. Alternatively, the satellite can send out 
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electromagnetic signals that are received at four locations (e.g., on the Earth surface) , and these four quantities can 
be used to compute the satellite position at emission time. 

In either case, the equations for navigation in a curved space-time can be formulated in a covariant, and invariant 
way (independent of coordinates) using the world function IJ. For example, consider a satellite with unknown 
coordinates = {to,x.o), z = 0, 1, 2, 3, which we want to locate with respect to four electromagnetic beacons having 
coordinates a;^ = (ia,Xa), a — 1,2,3,4. Assume that the satellite simultaneously receives the four signals at the 
space-time event with coordinates x^. We must solve the four equations given by 

n{xlxi) = 0, a = 1,2, 3, 4 (33) 

for the unknown satellite coordinates, = (to,Xo), in terms of the known emission event coordinates, xj, = (ta,Xa). 
These equations state that the emission and reception events are connected by null geodesies. In addition to Eq. H33|l . 
the appropriate causality conditions to > ta, for a = 1,2,3,4, must be added. The set of relations in Eq. H33|l are 
manifestly covariant and invariant due to the transformation properties of the world function under independent 
space-time coordinate transformations at point Pi and at P2. Equations l|33() neglect atmospheric effects, although 
these could be included in future work. 

From the definition of the world function, the intrinsic limitations of navigation in a curved space-time are evident: 
the world function Q{Pi,P2) must be a single- valued function of Pi and P2. In general, if two or more geodesies 
connect the points Pi and P2, then rj(Pi,P2) will not be single- valued and the set of equations in Eq. may 
have multiple solutions or no solutions. Such conjugate points Pi and P2 are known to occur in applications to 
planetary orbits and in optics |3lj |. However, when the points Pi and P2 are close together in space and in time and 
the curvature of space-time is small, we expect the world function to be single valued and the solution of Eq. H33|) 
to be unique. Therefore, navigation in curved space-time is limited by the possibility of determining a set of four 
unique null geodesies connecting four emission events to one reception event. In the vicinity of the Earth, there is no 
ambiguity, due to the weak gravitational field. However, in the case of strong gravitational fields, as may exist in the 
vicinity of a black hole, or when the (satellite) radio beacons are at large distances from the observer in a space-time 
of small curvature, navigation by radio beacons may not be possible in principle. In such cases, one may have to 
supplement radio navigation by inertial techniques; see for example the discussion by Sedov |6lj . 

V. PHYSICAL MEASUREMENTS 

In order to extract information from a quantum mechanical system, a measurement has to be performed. In a 
quantum mechanical theory, the measurement process plays a central role. The quantum mechanical measurement 
process is irnagined as an interaction that occurs between a classical, macroscopic apparatus and the quantum sys- 
tem ^3123, • Much effort has been expended on detailed investigations of the measurement process. 

In contrast, in the general theory of relativity, comparatively little attention has been devoted to the role of the 
measuring process. In part, this is due to the fact that really high-accuracy measurements were not carried out due to 
the associated technological difficulties. However, the presence of high accuracy clocks in space make high accuracy 
measurements a reality. 

A. Observations and Measurements 

In a real laboratory experiment, the process of taking data consists of recording information that is observed on 
physical meters, dials and gauges. For example, the reading on a volt meter can be observed at some point in time 
and space. The reading on the volt meter is an example of an observation. The observation occurs at some point 
in space and time (an event), marked by coordinates in some 4-dimensional system of space-time coordinates. The 
space-time coordinates of an observation may not be known to the person that is making the observation. 

Physically useful observations usually consist of coincidences of two (or more) things occurring at the same place 
at the same time. As an example of an observation, consider walking on a side walk, and noting the time on your 
watch when passing a crack in the sidewalk. The time noted is an observation of a coincidence: the position of my 
watch with the position of the crack in the sidewalk. This coincidence occurred at a space-time event, which has 
some coordinates. The observation is idealized as occurring at a given point in space-time, and each situation must 
be analyzed whether the "region of observation" satisfies, for the given application, the accuracy requirement of the 
observation as occurring at a point. 

In contrast to an observation, a measurement is an observation in which a comparison is made. For example, consider 
an a.c. electric current that is fed into a circuit where the phase of the incoming current is compared {measured) 
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rod to be measured 




O 

observer 

FIG. 1: The standardized meter stick is shown next to the rod to be measured. 

to the phase of a reference a.c. current. The measurement is the phase difference between these two currents. The 
measurement occurs at a specific space-time event, with definite coordinates. In the real world, the event is not a 
point, but often, to a sufficient degree of accuracy it can be modelled as occurring at one space-time point. Whether 
any given measurement can be regarded as occurring at one point in space-time depends on the required accuracy, 
and must be analyzed on a case by case basis. 

As another example of a measurement process, consider measuring a rod, by use of a standardized meter stick. 
Light from the ends of the rod comes to our eyes, along with light from the graduated scale on the standardized 
meter stick. The (simultaneous) event of light entering our eye from the left and right sides of the rod and meter stick 
constitute an observation, and because it is a comparison, it is also a measurement. See Fig. ^ Measurements are a 
subset of observations. 

Measurements are dimensionless ratios: the thing measured is compared to a standard. Furthermore, measurements 
are invariant quantities. In general relativity theory, the tetrad formalism treats measurements as invariant quantities. 

B. Tetrad Formalism 

Consider two observers that are making spectral measurements on light from the same star. Assume that the two 
observers are in relative motion, but that at the instant of measurement, they are located at the same point in space. 
At this point, the observer that is moving toward the star may measure predominantly blue light emitted from the 
star. On the other hand, the other observer that is travelling away from the star may measure predominantly red 
light. So two measurements at the same place at the same time lead to different results. (There are other types 
of measurements that may produce identical results for the two observers.) In the previous section, we stated that 
measurements are invariant quantities. In what sense then are measurements invariant? 

In connection with measurements, there are two types of transformations that must be considered. First, the global 
coordinates in the space, x', can be transformed to new coordinates, say using a transformation such as 

x^-^y^^ r{x^) (34) 

where /*() are a set of transformation functions. Since measurements are scalar quantities (see below), they are always 
invariant under the generalized coordinate transformations of the type in Eq. I|34|l . 

The second type of transformation that must be considered in connection with measurements is that two observers 
have different world lines and consequently different tetrad basis vectors onto which they project electromagnetic 
fields. The projection onto the tetrad is the measurement process. Real measurements are local quantities, and they 
can be compared when two observers are colocated at the same space-time point, see Figure |21 Transformations from 
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(1) 




l\o. = V ' 



FIG. 2: The world lines of two observers, A and fi, are shown with their respective tetrads, AJ^j and ^{q)- point M, 
the observers are momentarily colocated and they make a measurement of the same physical quantity. The results of their 
measurements are related by the transformation between tetrad basis vectors given in Eq. 13511 . 

the tetrad basis A^^^, a = 0, 1, 2, 3, of one observer to the tetrad basis of another observer, /ij^), can be contemplated: 

for a,/? =1,2, 3 (35) 

where H^^^^ is a 3x3 rotation matrix (a,/3 — 1,2,3) that relates the spatial tetrad basis vectors. Note that the 
matrix -ff (f j relates three 4- vectors A^^^ to three 4- vectors /i^^-j . (For each observer, the 0th components of the tetrad 
basis, A^Qj and ^^q^, are determined by their respective 4- velocity (see below), so these vectors do not enter the 
transformation . ) 

The key idea is that measurements are scalar quantities that are projections on the local basis vectors carried 
by each observer. Even though two observers coincide in time and space, their tetrad basis vectors are different: 
A^^^ for one observer and /i^^j for the other observer. Consequently, the two observers obtain different values of the 
measurement. Measurements are quantities that are projections on the observer's tetrad; they are of the form 

F^j>^)^ia) (36) 

for one observer and of the form 

F^, (37) 

for the other observer. Each of these quantities is a scalar, i.e., each is invariant under general coordinate transforma- 
tions given in Eq. 134|) . However, the measured quantities depend on each observer's tetrad and so the measurements 
are different. At any point, the observer's tetrad is determined by Fermi- Walker transport of the tetrad from some 
initial point, along the world line of each observer. 

The need for the tetrad formalism to relate experiment to theor y, a s well as the problem of measurable quantities 
in general relativity, is extensively discussed by Pirani '6^ , Synge |3l| , Soffel , Brumberg ,65] , and more recently 
within the context of metrology by Guinot ,66j . 

The tetrad formalism was initially investigated for the case of inertial observers that move on geodesies |63l l6^ IsM 
EM 1^1^' Many observers are terrestrially based, or are based on non-inertial platforms and the general theory for the 
case of non-inertial observers has been investigated by Synge [sj , who considered the case of non-rotating observers 
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movine; along a time-like world line, and by others [tIFIFI Fl FIFtIIt^ . who considered accelerated, rotating 
observers. Perhaps the most significant work for space-time navigation in the vicinity of the Earth is contained in 

Ref. [lEl Eliza III i3,E3 

As an illustration of the relativistic measurement process, consider an antenna that receives radio frequency elec- 
tromagnetic waves. The antenna converts an antisymmetric 4-dimensional tensor of second rank, Fij, into a scalar 
voltage reading on a meter. The meter may have a digital readout of the measurement. Consequently, the voltage 
is a scalar that does not transform under Lorentz (or generalized coordinate transformations). The voltage that is 
measured by a moving observer, V{t), is a function of the observer's proper time (since some starting point), t, and 
depends on the observer's tetrad defined on his world line. The voltage measurement process can be modelled as 

V{t) = F,,A'pi (38) 
= F.^XUKb^^M^"'^ (39) 

= ^^(a6)M('^^) (40) 

where Fij is the electromagnetic field in the space-time, and M'^ is the measurement tensor that depends on the 
world line of the observer. The scalar product of the tensors in Eq. H38|l reduces the electromagnetic field to a scalar 
quantity, V(t), which is the measured voltage. This voltage V{t) is invariant under coordinate transformations of 
the form in Eq. (PH. The projection of the electromagnetic field tensor, Fij, on the tetrad basis, A^^-j, yields a set 

of scalar numbers, -^(06) = ^ij A^^j A^^^^ . The quantity M^"^^ is an invariant matrix of numbers that characterizes the 
measuring apparatus (the observer's antenna). A simple model for an antenna is: 
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V '3 '772 -rni 



(41) 



where the six quantities, la and 777^, represent the sensitivity of the antenna to electric and magnetic fields, respectively. 
The 3-vector la is simply the vector effective length that characterizes the receiving antenna, see Eq. H27(l . 



1. Construction of the Tetrad: Fermi- Walker Transport 

The tetrad formalism for a given observer is constructed using the world line of the observer. Consider a time-like 
world line C of an observer given by coordinates x'^{u) with ui < u < U2, (see Fig. Along this world line, the 
observer has a four velocity = dx^ /ds and acceleration 

-^ = ^ = ^ + r>^.'= (42) 

OS ds 

where r*j. is the affine connection and Su^/5s indicates covariant differentiation along the world line x'(s). The 
normalization of the 4- velocity, u^Ui = — 1, provides the relation between the arc length, s = ct, where r is the proper 
time, and the parameter u. 

As the observer moves on the time-like world line C, he carries with him an ideal clock that keeps proper time r, 
and three gyroscopes. At some initial coordinate time x^, the observer is at point Pq at proper time r = s/c — 0. 
On his world line, the observer carries with him three orthonormal tetrad basis vectors, A^^^, where a = 1, 2, 3, labels 
the vectors and i = 0,1,2,3 labels the conipqnents of these vectors in some global system of coordinates. These 
vectors form the basis for his measurements [3lll63l. s ee Fig.|2| The orientation of each basis vector is held fixed with 
respect to each of the gyroscopes' axis of rotation |83|. The fourth basis vector is taken to be the observer 4- velocity, 
A^Q-) = 7x\ The four unit vectors A^^^, a = 0, 1, 2, 3, form the observer's tetrad, which is an orthonormal set of vectors 
at Po 

5y A(„) A^^) = r]ab (43) 

where the matrix rjf^^i,-^ is the Minkowski metric, see Appendix B. 

At a later time s — ct > 0, the observer is at a point P. The observer's orthonormal set of basis vectors are related 
to his tetrad basis at Pq by Fermi- Walker transport. Fermi- Walker transport preserves the lengths and relative angles 
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FIG. 3: The observer's world line C is shown with the initial tetrad basis vectors Aj^,) at s = at point Pq. The Fermi 
transported tetrad basis vectors at proper time s axe shown at point P. 



of the transported vectors. For an arbitrary vector with contravariant components its components at P are related 
to its components at Pq by the Fermi- Walker transport differential equations 

^ = W-'fj (44) 

where 

W'^ = u'w^ - w'u' (45) 

When we use Eq. 144|l to transport a vector that is orthogonal to the 4- velocity, u'/^ — 0, the second term in Eq. 145|l 
does not contribute. We refer to transport of such space- like basis vectors as Fermi transport, and W^^ W^^ — u^w^ . 
The space-like tetrad components satisfy A|^-|Mi = 0, for a = 1, 2, 3, and at any point P they are found by integrating 

the differential Eq. H44I) over the world line a;'(s), using the initial conditions in Eq. (|43|) on the tetrad at point P^, 
with W W''. 
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2. Fermi Coordinates 



Associated with each Fermi-transported tetrad basis, there is a set of Fermi coordinates, defined by the geometric 
construction shown in Figure |31 Every event P' in space-time has coordinates a;*' in the global coordinate system. 
According to the observer moving on a time-like world line, the same event has the Fermi coordinates X'^°'\ a = 
0, 1, 2, 3. The first Fermi coordinate, AT'^"' = s, is taken to be the proper time (in units of length) associated with the 
event P' . The proper time for P' is defined as the value of arc length s such that a space-like geodesic from point P 
passes through event P', where the tangent vector of this geodesic, /^*, is orthogonal to the observer 4- velocity at P: 

ti'u,\p^G (46) 

The orthogonality condition in Eq. H46|) is 

g.,^l\s)\\,^{s) = Q (47) 

and gives an implicit equation for s for a given point P' . This orthogonality condition gives the first Fermi coordinate 
of the point P' 

= s (48) 

The contravariant spatial Fermi coordinates, X^°'\ a = 1, 2, 3, are defined as 

= a^^'Af ) = g,, a{s) ^l\s) v^'^"^ Af^)(s) (49) 

where a is the measure along the space-like geodesic between P and P', gij is the metric and 77^*-'^ = //(yj = rjij is the 
invariant Minkowski matrix, see Appendix B. 



3. Metric in Fermi Coordinates 



All measurements made by a real observer are done locally, at the origin of Fermi coordinates. The measurements 
are projections on the tetrad of the observer, and the tetrad is only defined on the world line of the observer. Never- 
the-less, Fermi coordinates can be defined off the world line of the observer, and a corresponding metric for Fermi 
coordinates can be defined. The space-time interval in the Fermi coordinate system of the observer is 

ds^ = ~G^y) dX^'^ dX^j'^ (50) 

where G^ij) are the metric tensor components when the A"^*) are used as coordinates [3ll| . 

Despite the clear interpretation of measurement that the tetrad formalism offers, the analysis of experiments 
has seldom been done using the full tetrad formalism described above. In part, this is due to limitations of data 
accuracy and theorist patience to carry out the detailed computations. There are, however, a few explicit theoretical 
constructions of tetrads in the literature that address the issues discussed here [TII [tE [tM, Isa Ull • 



VI. REFERENCE FRAMES AND COORDINATE SYSTEMS 



The world function of space-time (see Section V) is a useful tool for describing the problem of navigation in space- 
time in a simple geometric way. The general equations are covariant, and invariant, and there is no need to specify a 
system of space-time coordinates. Navigation is carried out by making local measurements, in the comoving frame of 
the apparatus that does the measurement. 

However, in actual applications, a specific implementation of a system of space-time coordinates must be used. 
In most real-life applications today, such as the GPS, a fully relativistic 4-dimensional space-time coordinate system 
is not used. Instead, a system of three-dimensional coordinates plus a set of clocks are used. There are three 
common systems of three-dimensional coordinates that play a role in the navigation problem: Earth-centered inertial 
(ECI) coordinates. Earth-centered Earth-fixed (ECEF) coordinates, and topocentric coordinates. The origin of ECI 
coordinates is at the Earth's center of mass, and the orientation is determined with respect to distant objects-so the 
coordinates do not rotate with respect to these distant objects. However, the origin of ECI coordinates revolves around 
the sun along with the Earth, so these coordinates are better called quasi-inertial coordinates. ECEF coordinates have 
the same origin as ECI coordinates, but these coordinates rotate with the Earth, so a point that is stationary on the 
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Earth surface has a constant value for its ECEF coordinates. Topocentric coordinates have their origin on the Earth 
surface, with the x-axis pointing South, the y-axis pointing East, and the z-axis pointing radially away from Earth 
center (pointing up). Topocentric coordinates are used for making radar and other observations on the Earth surface. 
All three of these coordinate systems are three-dimensional. In experimental observations, these three dimensional 
coordinate systems are used together with clocks to record the coordinates of space-time events. Great care must be 
exercised when relativistic theories are used to analyze the data, because the observations were not recorded with 
respect to a true (relativistic) four-dimensional system of space-time coordinates. As an example of potential errors, 
see the next section. 



A. Gravitational Warping of Coordinates 



The presence of a gravitational field leads to a warping of the system of coordinates. More precisely, if we attempt 
to use Euclidean space and time coordinates, and neglect the effect of a gravitational field on the coordinate system, 
we will find that the definition of these coordinates is not precise. In other words, the definition of these Euclidean 
coordinates is ill-defined when a gravitational field is present. To demonstrate the magnitude this effect, consider the 
Schwarzschild metric 

- ds^ = g,jdx'dxi = (1 - "^ydt^ - - r2(sin2 Odcf^ + dO^) (51) 

^ r 

where the gravitational radius is given by Tg = 2GMJc^. For Earth, M =5.98 xlO^"* kg, so Tn = 0.88 cm. The 
physical length dl in a static space-time is given by 



IMJc 

m 



dp = dx'^dx^ = ^ + r2(sin2 ed<p^ + d9^) (52) 

r 

where the 3-dimensional spatial metric 70-/3 is related to the 4-diniensional metric gij by 

lal3 = (53) 

500 

The physical distance between two points at radial coordinates ri and r2 , and at the same coordinate angle and cj), 
is given by 



(r - rg) 



1 



+ rg log(v^ + - Tg) 



(54) 



Now imagine the first point is on the Earth's equator, so that ri = 6.378 x 10^ m is the Earth's equatorial radius, 
and the second point is at the altitude of a GPS satellite where r2 = 2.6561 x 10^ m. Since the Earth's gravitational 
radius is rg = 2GM/c? = 0.88 cm, the length I in Eq. can be expanded in the small parameter x = Tgjrx « 1, 
giving 



I = r2 — ri + - r„ loj 



ri 



(55) 



The first term on the right, r2 — ri, is the Euclidean length between the coordinate points at radial coordinate ri 
and at r2- The second term is the correction to the physical length I due to the presence of a gravitational field. 
The magnitude of this correction, for the ri and r2 values cited above, is ^Vg log ^ = 0.63 cm. This shows that the 
physical length I is longer by 0.63 cm than the difference between radial coordinate values in a Euclidean space that 
has zero gravitational field. The gravitational field has the effect of stretching the physical space between coordinate 
values. 

The implication of the effect is that a 4-dimensional relativistic frame of reference (coordinate system) must be 
implemented (rather than a system of 3-dimensional coordinates plus a time scale) if we want to have an unambiguous 
definition of a space-time coordinate system. In other words, a Euclidean system of 3-dimensional coordinates, plus a 
time scale, will have inherent error of ^Vg log ^ over a distance ^2 — ri, due to the warping of space-time due to the 
presence of the gravitational field. 
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In the above calculation, I have neglected the Earth's quadrupole moment J2 = 1.0826800 x 10~^ of the mass 
distribution. The inclusion of this quadrupole moment can be expected to modify the gravitational correction (second 
term) in Eq. H55(l by additional terms of order J2f'g, leading to an additional length correction on the order of 
J2rg log — , which is of the order of 10~^ mm over distances ^2 — ri = 2.0 x 10*" m. 

Another way of stating the effect of the gravitational warping of coordinates is to say that, if we used Euclidean 
geometry, the accuracy of clock synchronization is limited to (0.63cm)/c = 2.1 x 10^^^ s = 21 ps, even if perfect clocks 
and equipment are used. This value of 21 ps assumes the previous values of ri and r2. Reference to Table I shows 
that, for some applications, we must have time synchronization to better than 21 ps. 

The gravitational warping of coordinates is closely associated with (but distinct from) the slowing down of light in 
a gravitational field, also known as the Shapiro time delay effect. The average speed of light can be computed along 
a radial path from ri to r2. The path is a null geodesic given by ds^ = 0, which gives c^dt^ — dr"^ /{I — rg/rY . The 
coordinate time for light to traverse this path is 

1 dr 

dt^- (56) 



The average speed of light over this path is then given by 



I 



t2-ti 



1_ ^— log — 

2 r2 - ri ri 



(57) 



The second term in the above equation is the correction to the speed of light due to the presence of the Earth's mass, 
and for the previous values of ri and r2 has the small magnitude 

— = -i^^ log ^ = -3.13 X 10-1° (58) 
c 2 r2 — ri ri 

The negative sign in Eq. (|57|l means that light is has travelled slower than in vacuum in the absence of a gravitational 
field. This correction to the speed of light depends on the direction of light travel. 



VII. CLOCK SYNCHRONIZATION 



As discussed in the introduction, accurate clock synchronization is the backbone of applications such as high- 
accuracy navigation, communication, geolocation, and space-based interferometer systems. Also, as previously dis- 
cussed, clock synchronization cannot be divorced from the general problem of navigation in space-time (determining 
the position and time of an observer). The key point here, is that clock synchronization (divorced from the problem of 
determining spatial coordinates) is not a covariant concept, whereas navigation in space-time is a covariant concept, 
and consequently it can be formulated in covariant equations, such as was done in Subsection V-A, using the world 
function. See Eq. H33I) . However, if there exists a system of coordinates in which the clocks are all stationary at 
constant known spatial positions, then we can discuss their synchronization. However, this is generally not the case 
because satellites that carry clocks orbit the Earth. 

In the next Subsection B, we discuss synchronization of clocks that are stationary, at known spatial positions of 
an arbitrary system of coordinates, in the presence of a gravitational field. The system of coordinates need not be 
rigid and can be rotating. Section C discusses the peculiar way that GPS clocks are "synchronized" . Since the GPS 
satellites are moving, the clocks are (approximately) locally synchronized to a coordinate time in a certain metric, 
and hence the meaning of synchronization is different for GPS satellites than for stationary clocks. 

In pre-relativistic physics, clock synchronization was a straight-forward concept wherein two clocks were simply set 
to the same time. With the advent of highly accurate atomic clocks on board satellites, more accurate schemes for 
synchronizing remote clocks have been developed. Furthermore, the clocks to be synchronized are in relative motion 
as well as in different gravitational potentials. The accuracy of clocks has improved to the point that previously small 
unobservable effects must be accounted for by a comprehensive theory. A metric theory of gravity, such as general 
relativity, is the vehicle of choice for dealing with clock synchronization. Within the theory of relativity there are a 
number of clock synchronization schemes. However, one central concept-that of simultaneous event-is at the heart of 
the issue of clock synchronization. Simultaneity of two events is a definition, and not an absolute covarian concept, 
unless the two events are co-located at the same event in space-time. 
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A. Eddington Slow Clock Transport 

Clocks on maritime vessels of the 1700 and 1800's were used to navigate the oceans. An accurate clock was kept 
on the ship, and it was used to determine longitude of the vessel's current position [s^- The basic idea was to carry 
an accurate clock, and be in possession of the "correct time" . Using this time, and a view of the sky, the ship's 
longitude could be computed based on a theory of Earth rotation. This scheme of having the correct time by slowly 
transporting a clock is often called "Eddington slow clock transport" (84| . 

From the point of view of special relativity, where goo = —1 and gaf3 = Saf3, a, (3 = 1,2,3. (constant and uniform 
gravitational field) , a clock can be transported with arbitrarily small velocity and still maintain the "correct time" . 
From Eq. ((JJ, when the velocity of a clock is arbitrarily small, dx"^ /dx^ — > 0, a = 1, 2, 3, the relation between proper 
time T and coordinate time t = x^/c, reduces to At — Aa;'^/c. So a slow moving clock, which keeps proper time 
time T, can by definition be made to keep coordinate time . Such a clock can be slowly transported over large 
distances and can be used to synchronize other remote clocks. Unfortunately, Eddington slow clock transport is too 
restrictive for applications to clocks on satellites, because the speed of the satellites is not small, typically on the order 
oi v/c = dx"^ /dx^ ~ 10~^, and because the gravitational potential between different satellites varies significantly. See 
Section IX for a discussion of the effect of satellite motion and gravitational potential on proper time. 



B. Einstein Synchronization 



Another method of synchronizing stationary clocks is based on exchanging electromagnetic signals between clocks 
at known spatial locations. This type of synchronization is often called Einstein synchronization, and is based on a 
particular definition of simultaneity of (spatially separated) events. Perhaps the clearest discussion of the definition 
of simultaneity, within the context of general relativity, is given by L. D. Landau and E. M. Lifshitz We use 

their definition of simultaneity in the clock synchronization argument presented below. 

Consider two clocks at rest in some frame of reference (4-dimensional coordinates). Clock A has world line defined 
by constant spatial coordinates a;", and clock B has world line defined by spatial coordinates + c?x", a = 1, 2, 3, 
see Figure 0] Clock B is to be synchronized to clock A by exchange of electromagnetic signals. At event Pi, with 
coordinate time x^ + dx\, an electromagnetic signal is sent from clock B to clock A. Clock A receives the signal at 
event Pa, which has coordinates (a:°,a;"), and immediately reflects the signal back to clock P, where it is received 
at event P2 with coordinates (x° + dx^, x°' + dx"). This reflected signal (which is received by P at P2) contains the 
proper time reading, ta, of clock A at event Pa- 

The question then arises: what point P (coordinate time) on world line of clock B is simultaneous with event Pa on 
world line of clock Al Detailed consideration of this problem leads to the conclusion that there is no preferred way to 
deflne a point on world line B that is simultaneous with event Pa- Therefore, we arbitrarily take the midpoint between 
Pi and P2 as defining the point P that is simultaneous with Pa- For electromagnetic propagation, the space-time 
interval must vanish: 

-ds^ ^ g.jdx'dx^ (59) 

= gcf}dx"dx'^ + 2goadx^dx'^ + goo {dx^f = 

Solving for the coordinate time dx", we get two solutions, corresponding to the two directions of signal travel between 
x" and x" + dx°' : 



dx° = -— (goadx" - \/{goagoi3 - goo9ap)dx"dxP] 
500 V ^ / 

dx2 = — — (goadx" + J {goagop - googa/3)dx"dxl^] 
500 V ^ / 



(60) 
(61) 



Note that if dx2 > then dxf < 0. Therefore, the coordinate time of event P, which is simultaneous with event Pa, 
is given by (using the midpoint definiton of simultaneity): 

goadx"" 



900 



(62) 



Since events P and Pa are simultaneous (by definition), we consider clock A and B to be synchronized when their 
proper times are equal at these events: 

TB = TA (63) 
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FIG. 4: The world line of clock A and clock B are shown. Clock B is synchronized to clock A by electromagnetic signals. A 
signal from clock B is sent out at event Pi, the signal reflects off the face of clock A at event Pa, and returns to clock B at 
event P2. 



This corresponds to shifting the origin (epoch) of proper time. 

(2) 

The reflected signal from clock A arrives at B at the time .The proper time interval between event P and P2 is 
given in terms of coordinate time by the integral: 



(2) 1 



P2 



ds 



(64) 



Now, consider clocks A and B to be stationary, so that their spatial coordinates, and + dx", are constant in 
time in the chosen coordinate system. Then the proper time elapsed on clock A between point P and Pa is given by 



.(2) 



1 



TB 



V-.900 dx 







c JxO+:^{dx1+dx9^) 

^=\l {goagofi - googa0)dx°'dx'^ 

— Onn V 



V-ffoo 



(65) 
(66) 



where in the last line the quantities are evaluated at P, at coordinate time given by Eq. 162|) . Using the condition in 
Eq. Eq. can be written as 



(2) , ^ 



V^goo 



\J [gaagap - gm9ai3)dx°'dxP 



(67) 



Equation (|67|) gives the condition for clock B to be synchronized to clock A, in terms of quantities that are observed 
by clock B. Note that ta is the proper time on clock A, as observed by clock B. Equation H67|l gives the proper time 
that must be set on clock B at event P2 so that clock B is "synchronized" with clock A at the earlier event P. Note 
that the synchronization was done between clocks A and B that are separated by an infinitesimal spatial interval dx" . 
In practice, for small gravitational fields, this infinitesimal interval can correspond to large distances. For the case 
of a flat space-time with a Minkowski metric, goo = ~^,gap = 5a(3, Eq. H67|l gives r|j^^ = ta + I/c where the known 
spatial distance between clocks is I = gapdx^dx^ . 

Some comments are in order on the practical application of the synchronization condition in Eq. (|67|) . From 
Eq. H62|) . we know the space-time coordinates of point P = {x'^ — , -j-dx°'), and the synchronization condition 
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in Eq. (|67|l depends on the metric components gij . For high-accuracy clock synchronization, we must know the metric 
components gij in the coordinate system in which the clocks are at rest. For the accuracy needed in some applications 
(see TableOJ , the metric is not known to sufficient accuracy. An even more serious problem is that the spatial positions 
of clock A and B must be known a priori, i.e., the clock positions are not determined by the synchronization protocol, 
and the synchronization protocol depends on these positions. Finally, in the real world, satellites are in motion with 
respect to one another and there is no single (simple) coordinate system in which more than two satellites are at rest. 

In conclusion, it becomes evident that for most applications, the critical problem is not clock synchronization 
(correlating clock times), but that of navigation (correlation of clock positions and times). See Section V subsection A. 

Consequently, in real applications such as in the satellite system known as GPS, a different synchronization scheme 
is used. In the GPS scheme, satellite clocks keep (approximately) the coordinate time in the underlying ECI coordinate 
frame. We discuss GPS clock synchronization in the next section. 

In practice, Einstein clock synchronization is the basis of the applied technique known as two-way satellite time 
transfer (TWSTT) which in practice gives very accurate time synchronization between two points that are in common 
view of the same communication satellite 189, lOfl • 



C. GPS Clock Synchronization 

The Global Positioning System (GPS) is a U.S. Department of Defense Satellite System consisting of approximately 
24 satellites orbiting the Earth. The system consists of four satellites in each of six inclined orbital planes of 55°. 
Each satellite of the GPS carries atomic clocks, and sends pseudorandom code signals that are the basis of providing 
accurate time and navigation signals |i^^85, 86, 87] . The GPS uses a unique time synchronization scheme, wherein the 
clocks send time, encoded on chips of the pseudorandom code, and these chips (with time stamps on them) are received 
at Earth monitoring stations. The stations determine any clock corrections needed based on the broadcasts from the 
satellites. The GPS satellites move at a fast speed (approximetely w/c ~ 10~^) and the clocks suffer relativistic time 
dilation of approximately 7 fis per day. On the other hand, the satellites are at a high altitude, and the clocks run fast 
(as compared with coordinate time in the ECI frame metric), by about 45 /iS per day, due to the high gravitational 
potential in which they operate 0- The net effect is that, with respect to a clock on the Earth's surface, the clocks 
appear to run fast by approximately 38 fis per day. 

Each GPS satellite has an orbit that is approximately circular. Consequently, all the clocks on the satellites behave 
in roughly the same way-they run fast by 38 /is per day, or in terms of frequency, as observed from the surface of 
the Earth the oscillators run fast by 1 part in 4.4xlO~^°, see Section IX, subsection E. If something were not done to 
counter this effect, all GPS satellite clocks would appear to run fast. Consequently, a "factory offset" is applied to the 
frequency of each clock oscillator (in software) in the amount -4.4x 10"^". The clocks on board the GPS satellites then 
keep (approximately) coordinate time in the ECI frame, and the clocks appear to (approximately) keep correct time 
as seen from the Earth's surface. This "factory offset" does not compensate for the different orbital eccentricities of 
the satellites, since for each satellite. The factory also does not compensate for the relative velocity between satellites 
and different moving GPS receivers, see Ref. f3| for a detailed discussion of the role of the Doppler effect in GPS 
measurements . 

The "factory offset" applied to the GPS satellite oscillators can be understood as another scheme for clock synchro- 
nization. There exists a coordinate time t — x'^/c, in some metric around the Earth, see Eq. I|()8|) . The coordinate 
time is a global quantity, as compared to the proper time, which depends on the given world line. By definition, clocks 
on the Earth's geoid keep coordinate time in the ECEF frame. In order to synchronize the GPS satellite clocks to 
coordinate time in the ECI frame, the "factory offset" is applied. With this offset, the GPS clocks keep approximate 
coordinate time in the ECI frame of reference. As each satellite moves through space, we can imagine that it passes 
"hypothetical" coordinate clocks (that are stationary in the ECI frame and keep ECI coordinate time). The "factory 
offset" has the effect that a GPS satellite clock instantaneously agrees (closely) with each ECI frame coordinate clock 
that it passes, see Figure^ So the GPS satellite clocks are (approximately ) synchronized to coordinate clocks in the 
underlying ECI frame. Actually, this is an unusual idea-for (GPS) clocks in one frame (comoving frame of a GPS 
satellite) to keep coordinate time in another (ECI) frame. As mentioned above, in the case of GPS, this synchroniza- 
tion is approximate, because only a constant rate offset is applied before satellite launch, which cannot compensate 
for effects, such as the eccentricity of each satellite orbit (called esin£^ effect), the Earth's quadrupole moment (or J2 
effect), the effects of other forces (such as solar pressure and atmospheric drag) on each satellite, and imperfect atomic 
clocks. The satellite ground tracking system for the GPS monitors satellite clocks and attempts to compensate for 
these effects by calculating satellite clock corrections and uploading these corrections to the satellites so that users of 
the GPS may apply them to get accurate space-time navigation. 
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FIG. 5: The motion of a GPS satellite in 3-dimensional space is shown as it passes (stationary) coordinate clocks in the ECI 
frame. At any instant, the GPS satellite clock agrees with time on the underlying ECI frame coordinate clocks. 

D. Quantum Synchronization 

During the past several years, alternative schemes for clock synchronization [21I |2^ 12^ 123. 125l j2^j2^ I2M I29I l30l| 
have been proposed based on quantum information theory and entanglement of quantum states [l9l |20| . Quantum 
effects may be exploited for clock synchronization since entangled (or correlated) photon pairs are found to interfere 
destructively at a beam splitter |2lj |. Such photon pairs are believed to be created at a single space-time event, 
within 100 fs of each other '2l'|. Exploiting this effect, two relatively simple synchronization schemes have been very 
recently proposed based on the interference of entangled photon pairs states, which can be created by parametric 
downconversion in a crystal that lacks a center of inversion symmetry . One of these schemes is capable of clock 
synchronization in free space '2^, while the other relies on there being a difference of group velocities for each of 
the photons in the entangled pair in an optical medium, such as exists in an optical fiber 30]. Ve ry r ecently, an 
experimental demonstration of quantum clock synchronization has been carried out in the laboratory jSQ] . 

The quantum mechanical clock synchronization proposals do not included the basic relativistic effects present when 
clocks on board satellites are to be synchronized: the fast relative motion of clocks and the variation in gravitational 
potential between satellites in various orbital regimes, such as low Earth orbit (LEO), geosynchronous orbit (GEO), 
and highly elliptical orbit (HEO). In the next section, we describe the resulting syntonization problem that arises for 
satellites at differing gravitational potentials. 

VIII. THE SYNTONIZATION PROBLEM 

Satellites in multiple orbital regimes may have their clocks synchronized by means of exchange of optical signals, 
as described in the previous section. This means that at one coordinate time, all clocks can be made to read the 
same value (for example, the same coordinate time). Satellite applications require that the synchronization between 
members of the satellite ensemble be maintained for some time period, or alternatively, the clock differences must be 
known at a given elapsed time from the synchronization epoch. It is well known that the proper time, and consequently 
the hardware time (see Section III), will run at a different rate on each satellite with respect to coordinate time of some 
given metric. This difference in rate of proper time is due to the motion of the clock (time dilation) and gravitational 
potential effects (sometimes called the gravitational red shift effect). 

Since the hardware time on all real clocks is affected by their motion and position in the space (local gravitational 
potential), the only reasonable measure of time is coordinate time, as used in a metric theory of gravity, such as general 
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relativity. Coordinate time is a mathematical construct that is global, which means that it is the same everywhere in 
the space. In distinction, proper time depends on the history (world line) of the clock. In this Section, we compute 
the difference between proper time on board a satellite, and coordinate time, for satellites in various orbits about the 
Earth. 

Since the clock on board a satellite is located at a different spatial position than a reference clock on the Earth's 
surface, these clocks can be compared by exchange of electromagnetic signals. The relation between these two clocks 
must be established, such that their relative motion and their respective gravitational potentials are taken into account. 
At any instant in time, the relative velocity of the satellite with respect to the ground leads to the satellite clock 
running slower than the ground clock, due to relativistic time dilation. The higher gravitational potential at the 
satellite leads to the satellite clock running faster than the ground clock. Generally speaking, clocks in low-orbiting 
satellites run slow due to the predominant time dilation effect, due to high orbital velocity and small altitude above 
the Earth's surface. On the other hand, clocks in high-orbit satellites generally run faster than ground clocks because 
the gravitational potential effect is predominant. 

While a satellite clock can be compared to an Earth-bound reference clock, there is nothing special about the 
Earth-bound clock as a reference clock. In fact, from the point of view of an ECI coordinate system, the Earth bound 
clock moves in a circle, just like the satellite clock. Furthermore, the actual reading on an ideal clock depends on its 
world line, i.e., its past history of velocity and gravitational potential. All ideal (and hardware) clocks suffer from this 
complication. Therefore, the comparison of all clocks must be made to a standard time or "clock" that is global and 
does not depend on its history. Such a quantity is the coordinate time associated with some metric that describes 
the space-time in the vicinity of the Earth. Coordinate time is a mathematical construct that is the same in all of 
space-time. The complication is that coordinate time (and generally each spatial coordinate) is arbitrary in a metric 
theory such as general relativity, and depends on an arbitrary choice of a space-time coordinates. In the next section, 
we will choose a space-time metric that has the desirable property that coordinate time corresponds (approximately) 
to the proper time kept by an ideal clock on the Earth's geoid (surface of equal geopotential). 



The definition of coordinate time comes from choosing a specific metric for the space-time. In general relativity, the 
coordinates x*, i — 0,1, 2, 3, are mathematical entities that are never observed. The metric of space-time depends on 
these coordinates, and takes into account the warping of space-time due to the presence of the gravitational field. For 
the same physical space-time, we can choose different coordinates. Of course, all observations are independent of the 
coordinates, and so the choice of coordinates is arbitrary. However, it is convenient to choose space-time coordinates 
so that coordinate time has some physical meaning. Perhaps the most reasonable choice is to take the metric for 
space-time of the form ^88j 



where x'^ , , x'^ , are geocentric equatorial coordinates, where x^ coincides with the Earth's axis of rotation and 
increasing positive values point to North. The Earth is taken to be an oblate spheroid with potential given by (9]| 



where, — {x^)'^ + (x^)^ + (x^)'^ and 6 is the polar angle measured from the a;'^ axis. In Eq. (|69|l . G is Newton's 
gravitational constant, M is the mass of the earth, P2{x) — (3a;^ — l)/2 is the second Legendre polynomial, 
is the Earth's equatorial radius, and J2 is the Earth's quadrupole moment, whose value is approximately J2 = 
1.0826800 X 10^"^. The metric in Eq. is the solution [32| of the linearized version of Einstein Eqs. 

The coordinate time can be given a simple interpretation. Transform the metric in Eq. (|68|l to rotating ECEF 
coordinates y', using the transformation 



A. Choice of Metric in Vicinity of the Earth 



(1 + ^F)(dT")2 + (1 - If) [{dx')^ + {dx^f + [dx^f] 



(68) 




(69) 




(70) 
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Note that the coordinate time in the rotating frame, = a;°. In these ECEF rotating coordinates, the metric is given 
by 



ds^ = - 



9V 9V 



[dff 



- y'dy')df + {dy^f + [dy^f + [dx^f 



(71) 



From Eq. (|71|l . we see that stationary clocks in the ECEF coordinates (clocks which satisfy dy"' = 0) that are at the 
same value of geopotential 4>, where 



^V~'-nmy'f + {y^f) 



(72) 



have the same rate of proper time with respect to coordinate time y". The term geopotential is used because 
(j) takes into account the effect of angular velocity SI of Earth rotation. We have neglected the small cross-term 

Using the observation that clocks at a constant value of geopotential run at the same rate, it is advantageous to 
define the new coordinate time t 



1 



where (po is the geopotential on the Earth's equator: 

GM 



(73) 



(74) 



The dimensionless magnitude of this term is (jio/c^ — —6.96928 x 10 Using the transformation in Eq. (O, the 
metric becomes 



ds' 



(0 - M 



o^2 



1 - 



2V 



2%'dy'-yW)dy° 
c 



2V 



+ {l-—)[{dy'f + (dyy + idx'f 



(75) 



Eq. (|75|l gives the space-time metric in ECEF rotating coordinates y*. Note that an ideal clock that is stationary in 
EFEC coordinates (with dy" — 0), has proper time 



dr = ds/c 



-,1/2 



dy' 



(76) 



When this clock is located on the geoid, then — (j>o, and dr = dy'^/c, so this ideal clock keeps coordinate time. 
Hence a good hardware clock on the geoid can be used as a standard to keep coordinate time, x'^, in the rotating 
ECEF coordinates. Note that by Ea. H7U|l the coordinate time in rotating ECEF coordinates is the same as coordinate 
time in ECI coordinates, so this same clock keeps coordinate time, x^, in the ECI frame coordinates a;\ 

Using the coordinate time transformation in Eq. (|73|l . the metric in the ECI coordinates, given in Eq. H68|l . becomes 



ds = gijdx^dx-' = — 



[dxy 



[{dx^ + {dx^ + [dx^] 



(77) 



Equation 1)77(1 gives the metric in ECI coordinates. The coordinate time that enters into the metric, x", is the time 
kept by ideal clocks on the geoid. This result was the singular goal of the time transformation given in Eq. ((73(1 . 

Note however, that in the ECI metric in Eq. 1(77(1 . the proper time interval ds on a stationary clock in ECI coordinates 
(with dx" = 0), is not equal to coordinate time interval dx'^ because in general V ^ (po- 

The ECI coordinate metric, given in Eq. ((77(1 . will be used below as the basis for comparing satellite clocks in 
different orbits. The time kept by ideal satellite clocks (proper time) will be compared to the global coordinate time 
x", which is the time kept by ideal clocks on the geoid. 
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B. Integration of Geodesic Equations 

In order to compare the proper time on a satellite clock to the coordinate time in ECI metric given in Eq. H77() , the 
satellite world line must be known. The satellites mass is negligible compared to that of the Earth, so the satellite 
is essentially a freely falling test particle that follows a geodesic in 4-dimensional space-time. The geodesic equations 
are given by 

''^ V,j2Li^,, (78) 



ds ds 

where the coordinates are given by = (a;", x°'), a = 1, 2, 3. The four velocity can be written as 

dx^js) dx\ dx'- dx^ , 

where v"^ — {l,dx°'/dx'^) — {l,dx"/d{ct)). The 4-velocity is normalized, 

- 1 = f ^ j [-a + pS^^v'^v^] (80) 

where a and /3 are functions of s. 

The four geodesic differential Eq. 1)78(1 must be supplemented by initial conditions: initial 4-velocity v*(s — 0) and 
an initial 4-dimensional position x'(s = 0). These initial conditions are taken from the classical Newtonian orbital 
mechanics in 3-dimensional Euclidean space. 

Now introduce Cartesian perifocal coordinates (z°, z^, z^, z"^), defined so the ori gin c oincides with one focus and the 
-|-z;"'^-axis contains the periapsis, and the z^-z^ plane contains the elliptical orbit '92'|. In these perifocal coordinates 
z', the Cartesian coordinate position r and velocity v are given by 



r cos i^P -1- r sin z^Q (81) 
[- sin z/P -f (e -f cos v) Q] (82) 



GM 



a(l-e2) 

where v is the true anomaly, a is the semimajor axis of the orbit, and e is the orbital eccentricity. In Eq. (|81() - (|82|l . 
P and Q are three-dimensional orthogonal unit basis vectors along the z^- and z^-axis. In perifocal coordinates, the 
radial coordinate of the orbit is given by 

. = ^(i^ (83) 
1 + e cos V 

For the initial conditions, take = 0, so the Cartesian position and velocity components are given by 

r(0) = (zi(0),z2(0),z3(0))-(a(l-e),0,0) (84) 
' dz^ dz^ dz^^ 



v(0) 



Note that in Eq. (|85|l I have taken the origin of coordinate time z'^ = s — 0, when the satellite is at the periapsis. 
This is an arbitrary choice made for convenience, and does not affect the results computed over one orbital period. 
The transformation between the perifocal coordinates, z*, and the geocentric equatorial coordinates, x", is given 

by 

= zO (86) 
a,/3=l,2,3 (87) 



c — d ^ z^ , 



where the orthogonal transformation matrix is given by 

(cos(r2) cos(tj) — sinfi sinw cosi —cos $7 sin a; — sinfi cosw cos« sinf2 sini 

sin f2 cos Lu + cos sin uj cos i — sin f2 sin uj + cos ft cos uj cos i — cos 51 sin i 

sinw sini cosw sini cosi 
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The 4-dimensional initial conditions are taken as 

x°{0) = (89) 

x°'{0) = d"^ z'3(0) (90) 

w°(0) = 1 (91) 

-"(0) - d-, ^ (92) 
Using the metric in Eq. H77|l . and the normalization of the 4- velocity 

- 1 = 5,,vV^ (93) 



I find that 

ds 

where 



/3(5,^t;''i;^]-^/' (94) 



a = l + 4(V^-'/'o) = l + ^ (95) 

/3 = 1 - 4-^^ = 1 - (5 (96) 
c/ 

and both a and [3 are functions of s through the coordinates x°'{s). 
The 4- velocity can then be expressed as 

V{s) = [a ~ pS.^v'^v^]-^'^ (!,«") (97) 
So the initial condition on the 4-velocity become 

v^O) = {[a - (35,^v''v^y^'^}s=o (1, f "(0)) (98) 

where w"(0) is given by Eq. (|92|l . In all formulas, Greek indices take values a — 1,2, 3, and Latin indices i = 0, 1, 2, 3. 

Using Eq. (|80|l and H79|l . the rate of change of satellite proper time with coordinate time can be written in terms 
of the 4-velocity vector components as 



ds 
d^ 



1 + pS^^v^v^ 



-,1/2 

, K,7= 1,2,3 (99) 



Equations H78|) can be integrated using the initial conditions given in Eq. (|89|I - H92() . The satellite orbital elements 
enter the differential equations through Eq. l|5H) - (|S3|) and Eq. (|S5|l . 

The goal is to compute the difference between satellite proper time and coordinate time. We also compute the 
frequency shift (gravitational plus Doppler shift) as observed on Earth, from a satellite transmitting at a known 
frequency. The absolute orientation of the satellite orbit is not of interest here. Therefore, we set the longitude of the 
ascending node H. — tt/2 and the argument of periapsis u; = 3tt/2. Setting these two orbital parameters in this way 
puts the apogee on the negative x— axis, on the z > half-space. The perigee is on the negative x— axis, at z < 0. 
See Figure El 



C. Proper Time Minus Coordinate Time for Various Orbits 

The difference between proper time r kept by the satellite clock, and coordinate time At ~ Ax^ /c, kept by clocks 
on the geoid (in metric given by Eq. ifzTjl ) is given by 

1 dx° 

At{s) - At{s) = - / (1 - -—)ds (100) 

c J,^ ds 

where dx^ /ds is given by Eq. (|99|l . Equation (|100|l gives the time difference, At(s) — At{s), at a given value of proper 
time measured from some epoch, or starting time. We choose the arbitrary starting time si — ct — ct = x'^ = to be 
when the satellite is at periapsis. 



27 




FIG. 6: The orbit orientation is shown in a real-space 3-dimensionaI view at time x = t = 0, with perigee placed at z < 
and at a maximum value of the x-axis of the orbit. (The perigee is in the lower right side of the figure.) At time x'^ = t — 0, 
the receiver (shown as cube) is on the equator on the surface of the Earth at (x, y, z) = [R, 0, 0), where R is Earth's equatorial 
radius. 



The gravitational field of the Earth is such that the satellite orbit does not close on itself. In the integration of 
the differential equations for this model, we have arbitrarily chosen s = si = as the point at which a satellite is at 
periapsis. We take the orbital period to be the value of s = Sp at which the satellite is at its closest approach, as 
defined by the minimum of the Euclidean distance squared cP{s): 



(101) 



For satellites orbiting Earth, we have the small the quantities v « 1, 6 « \ and (5 << 1, and the integrand in 
Eq. HIOOI) can be approximated by 



l-^^if.~-.^)-i 
ds 2 V / 8 



3<52 -v"- AuH - 2vH 



0(5) 



(102) 



where we have dropped terms that we call 0(5). Terms of 0(1) ~ 10 ^, and so ^ 10 ^ ^ 0(1), (5 ~ 10 ^" ~ 0(2), 
(5 ^ 10"^" ~ 0(2), for typical Earth orbiting satellites. Terms in the integrand in Eq. (|100|l have been dropped that 
are of order 10~^^. For a satellite with an orbital period that is as large as one day (86400 sec), the error in dropping 
these terms is on the order of lO^^*' sec. So our results, within this model, are accurate to lO^^*' sec per day. 
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TABLE II: Numerical value of the parameters used in the computer program. 



c (m/s) GM (mVscc^) Jz 


R (m) 


Earth Rotation Rate Q, (radian/s) 


2.99792458x10** 3.986005x10" 1.0826800x10"^ 


6.378137x10" 


7.2921151467x10"" 


TABLE III: For satellites with given orbital parameters, semimajor 
quadrupole moment J2, the computed values of proper time minus 
and per day. The value per day is the average of this difference over 


axis a, eccentricity e, 
coordinate time, At - 
one Earth solar day. 


inclination i, and value of Earth's 
- At, are shown per orbital period. 



satellite 


semimajor axis 
(meters) 


eccentricity 


inclination 
(degrees) 


J2 




At period 
(minutes) 


At - At 
/^s per period 


At - At 
/is per day 


LEO" 


7.3635 xlO'^ 


0.00292 


82.9 


1.0826800x10" 


-3 


105.12 


-1.290509 


-17.678433 


LEO' 


7.3635x10^ 


0.00292 


82.9 


0.0 




104.81 


-1.301039 


-17.875853 


GEO" 


4.2164174 xlO^ 


0.0 


0.0 


1.0826800x10" 


-3 


1435.96 


46.4512489 


46.5818860 


GEO'* 


4.2164174 xlO^ 


0.0 


0.0 


0.0 




1436.0 


46.4230537 


46.5501514 


HEO" 


2.70365x10^ 


0.747194 


62.8 


1.0826800x10" 


-3 


743.08 


20.1582623 


39.0644760 


HEO^ 


2.70365x10^ 


0.747194 


62.8 


0.0 




737.37 


19.9308525 


38.9226991 


GPS" 


2.66965x10^ 


0.0017418 


55.03 


1.0826800x10" 


-3 


723.573310 


19.438916 


38.6858366 


GPS'' 


2.66965x10^ 


0.0017418 


55.03 


0.0 




723.504421 


19.420036 


38.6519441 



"COSMOS 2366 Russian military spacecraft 
''Same as above LEO but taking J2 = 0. 
■^Nominal values for a GEO satellite. 
■^Same as GEO above, but with J2 = 0. 
^Highly Elliptical Orbit: 1998-054 Molniya 1-91 
■'^Same as above HEO but taking J2 = 0. 
finternational Designator 1996-041 
''Same as GPS above, but with J2=0 

D. Proper Time Minus Coordinate Time: Numerical Results for Various Orbits 

A computer program was written in the Mathematica programming language to solve the general relativistic geodesic 
equations of motion given in Eq. (|78|l . subject to the initial conditions on the 4- velocity in Eq. (|98|l . The purpose of 
this program is to illustrate the relativistic effects on satellites in various orbital regimes, and not to predict accurate 
orbits that include all satellite perturbations. Consequently, the effects of atmospheric drag and solar pressure have 
been neglected, but in the future these effects could be included. The computer program takes as input the satellite 
orbit semimajor axis a, eccentricity e, and inclination i, as well as the gravitational constant times the Earth mass, 
CM, and the value of the earth's quadrupole moment, J2. The values of these parameters are shown in Table Hll 

As mentioned above, the orbit is arbitrarily taken to have the longitude of the ascending node 17 = 7r/2 and the 
argument of periapsis uj = 37r/2. Choosing these parameters in this way puts the apogee on the negative x— axis, on 
the z > half-plane, and the perigee on the negative a;— axis, at z < 0. The computer program solves the geodesic 
equations of motion numerically for the (equatorial geocentric) coordinates 2:*(s), z = 0, 1,2,3, as a function of the 
proper time s along the orbit. The equations of motion arc integrated from s = at perigee, to s = Sp at the closest 
approach to the initial position (approximately one orbital period) which defines the orbital period. At perigee, both 
the coordinate time and proper time are arbitrarily taken to be equal, and are set to zero. As the satellite progresses in 
its orbit beyond perigee, proper time deviates from coordinate time, due to relativistic time dilation and gravitational 
potential effects. 

In order to illustrate the numerical values of the effects discussed above, we use a representative satellite from each 
of four orbital regimes: low-Earth-orbit (LEO), geostationary orbit (GEO), Global Positioning System (GPS) satellite 
orbits, and highly-elliptical orbit (HEO). The numerical values used in the computer program are shown in Tabic ITTT 
These values were obtained from CelesTrak WWW on the the World Wide Web from web site: http:/ /celcstrak.com 
We compute the difference between elapsed proper time and coordinate time. At — At, where At is the elapsed proper 
time, and At is the elapsed coordinate time, at which the satellite is at perigee. We arbitrarily choose the coordinate 
time and proper time to be zero when the satellite is at perigee. Then, up to an additive constant, the coordinate 
time, t = x^/c, where appears in the metric in Eq. H77|l . is equal to the proper time kept by a reference clock on 
the geoid, such as the Master Clock at the U.S. Naval Observatory in Washington, D.C. 

For each satellite. Table Ulll shows the the computed values of the proper time minus coordinate time difference. 
At — At, from perigee to (approximate) perigee, accumulated over one orbital period, as determined by solving the 
geodesic equations of motion, and evaluating Eq. (|100() . The periods of the satellites vary widely, so to compare 
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FIG. 7: The proper time minus the coordinate time, Ar(s) — At(s) is shown for the LEO, GEO, HEO (Molniya) and GPS 
satellites vs. s/sp = r/rp, the fraction of proper time period along the orbit. On the abscissa, to the accuracy of the plot, 
proper time difference is approximately equal to coordinate time difference, Ar ~ At, where t is coordinate time measured in 
seconds. 



the magnitude of the effects, the last column in Table IIIII shows the difference Ar — At integrated over one solar 
day. This value is simply taken to be the average of Ar — Ai over one satellite period, multiplied by the ratio (24 
hour) /(satellite period). The calculations in the Table ITTll are done for a finite value of Earth quadrupole J2, and for 
J2 = 0. Comparing the time difference Ar — At for J2 = and for finite J2 gives an estimate of the size of the effect 
of the Earth's quadrupole on proper time kept by a satellite clock. The coupling to the quadrupole terms is larger for 
lower altitudes (since the terms proportional to J2 vary as l/r'^) and depends on the inclination of the orbit. 

For the low-Earth-Orbit (LEO) satellite, the apogee is at altitude of 1007 km and perigee at 964 km. The difference 
between proper time and coordinate time, Ar — At = —17.678433 ^s per day, which means that less proper time r 
(aboard the satellite) elapses over one solar day than coordinate time t on the geoid . This is due to the predominant 
effect of time dilation resulting from the high orbital speed of a LEO satellite. Also shown in Table ITTll is the effect of 
the non-zero value of the Earth's quadrupole moment, J2. For the LEO satellite, the finite value of J2 has the effect 
of changing the difference Ar — At by 197.42 ns. 

For high altitude GEO satellite, the difference between proper time and coordinate time, Ar — At, is of opposite 
sign to that of a LEO satellite. For a GEO satellite, the proper time elapses faster than coordinate time by 46.5818860 
jjLS per day. The gravitational shift (causing proper time to run fast) dominates the time dilation (which causes proper 
time to run slow), due to the satellites high altitude and relative slower speed. Also, the GEO satellite is less influenced 
by the finite size of Earth's quadrupole moment J2, which is only 31.7346 ns per day, as compared to 197.42 ns for 
the LEO satellite. Furthermore, for a GEO satellite, the contribution to Ar — At from the J2 quadrupole effect is 
opposite to that for a LEO satellite. 

For the (high-altitude) GPS satellite, apogee is at 20365 km altitude and perigee is at 20272 km altitude. The 
predominant effect is due to change of gravitational potential, which leads to more elapsed proper time than coordinate 
time, by 38.6858366 /its per day. Since the altitude of GPS satellites is lower than GEO satellites, the increase of 
proper time is smaller. For GPS, the effect of nonzero value of J2 is also similar to the GEO satellite, but has the 
larger value 33.8925 ns, due to the lower altitude which leads to a greater influence of J2 

The LEO, GEO, and GPS satellites have nearly circular orbits. In contrast, the HEO satellite has an eccentricity 
of 0.747194, and consequently a perigee of 6835 km (457 km altitude) and apogee 47238 km (40,860 km altitude). 
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FIG. 8: The proper time minus the coordinate time difference, Ar(s) — At{s), is shown for the LEO, GEO, HEO (Molniya) 
and GPS sateUites vs. coordinate time At in units of hours. Coordinate time At is equivalent to time kept by an ideal reference 
clock on the Earth's geoid. 



This places the HEO satellite both lower and higher than the nearly circular orbit LEO satellite during different 
portions of its orbit. The altitude of the HEO satellite at perigee and apogee is also lower and higher than the GPS 
satellite altitude. Consequently, the HEO satellite has proper time running more slowly than coordinate time when 
it is at low altitude (like the LEO satellite) and it has proper time running faster than coordinate time when its is 
at high ahitude (like GEO and GPS sateUites). The net effect for the HEO sateUite is that 39.0644760 /is per day 
more proper time has elapsed than coordinate time. This size of this net effect is similar to the GPS satellite because 
the HEO satellite spends little time at the low altitude portion of its orbit. The effect of J2 for the HEO satellite is 
141.777 ns, which is 4 times the same effect for the GEO satellite, and 0.7 times this effect for the LEO satellite. 

The results described above are the integrated value of proper time minus coordinate time, over one complete 
orbital period, ft is useful to plot the value of proper time minus coordinate time as a function of the fraction of 
the orbital period covered, s/sp, see Figure [7| This figure shows that for each satellite, the proper time At — s/c 
diverges from the coordinate time At — Ax^/c, at a different rate, along its orbit. Note that time is scaled to the 
orbital period along each orbit, so that s/sp = 1 occurs at the respective period of each orbit. The plot in Figure[7| 
shows a geometrical comparison, however, it does not permit a comparison of the effect in time. 

For some practical purposes, it is preferable to look at the proper time minus coordinate time difference. At — At 
as a function of coordinate time At — Ax^ / c. In this way, the effects are shown as they would appear in real time, say 
with respect to the USNO Master Clock. Figure shows At- At vs. At for the LEO, GEO, HEO (Molniya) and GPS 
satellites. For each satellite, on the scale of this graph, the plots are almost linear functions. The obvious exception is 
the HEO (Molniya) orbit satellite, where, due to its high eccentricity, the time difference is highly nonlinear. Actually, 
even for the LEO, GEO, and GPS satellites, which have nearly circular orbits, the rate of change of Ar(s) — At{s) 
with time t is non-constant, due to the interaction of the orbital eccentricity, inclination, finite value of J27 and general 
relativistic precession effects. These effects are not visible on this graph, but can be significant, depending on the 
desired level of time synchronization that is required (see Tabled and the time interval over which the synchronization 
must be maintained (stability). The nonlinear variation of proper time with satellite time can best be seen on a plot 
of the frequency shift, see the next Subsection. 

Since the rate at which proper time diverges from coordinate time is nearly constant for the LEO, GEO, and 
GPS satellites, these values could be used to compensate for the drift of proper time away from coordinate time, by 
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FIG. 9: The world line of the satellite (transmitter) and observer (receiver) are shown schematically. During satellite proper 
time dss, the satellite clock sends A'^ uniformly-spaced pulses, which are received during observer proper time dso- A typical 
reception event O is connected to an emission event S by a null geodesic. 



applying a "factory offset" to tlie satellite clock frequency, as has been done in the GPS. This could even be done after 
the orbit is established and the compensation could even be changed periodically. However, inspection of the time 
synchronization that is required for various applications, see for example Tabled shows that for many applications 
this compensation is not feasible, due to the non-constant rates (for reasons mentioned above), as well as due to other 
forces (such as solar pressure and atmospheric drag) that perturb the satellite orbit and lead to proper time changes 
that cannot be modelled sufficiently accurately. Instead, it is likely that for the most accurate applications of time 
synchronization, the synchronization will have to be redone periodically, according to the required stability (accuracy 
required over given time interval). 



E. Observed Doppler and Gravitational Frequency Shift 



When a satellite clock is manufactured, it has a certain prescribed hardware clock rate compared to coordinate 
time, in some system of coordinates. Alternatively, we can say that the clock, or oscillator in the clock, is calibrated 
to a given frequency lOs when it is compared to a local frequency standard. The word local here means that during 
the calibration procedure, the satellite clock is siting at rest with respect to the frequency standard and that these 
two devices are at the same gravitational potential (they are co- located and at rest relative to each other). We then 
refer to the frequency ujs as the proper frequency of the satellite oscillator or clock. 

When the satellite is launched, its clock is assumed to operate normally. In fact, if the frequency standard were 
placed aboard the satellite, and a comparison was made between the satellite clock and the frequency standard, we 
would expect to obtain ujs for the satellite oscillator frequency, to within the calibration accuracy. However, if the 
frequency standard is located on the Earth's geoid, and the satellite clock sends out electromagnetic pulses, then the 
frequency of the pulses observed on the geoid (with respect to the frequency standard on Earth) will be Wq, which 
is different than the proper frequency of the satellite clock, Wg. This frequency shift is due to the relative motion of 
satellite and receiver (Doppler effect) and also due to the difference in gravitational potential between satellite and 
receiver. In what follows, we compute this effect. 

Consider a satellite orbiting Earth on world line xl{ss), i — 0, 1,2,3, where the functions xl{) specify the world 
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line of satellite s, and Ss is the parameter related to proper time along the world line, t = Ss/c. See Figure Now 
consider an observer on a different world line, a;* (sq). The satellite broadcasts a periodic electromagnetic signal, and 
within the geometric optics approximation, we say that these signals travel along null geodesies, so the reception event 
O and transmission event S are connected by a null geodesic. In terms of the world function we have 



n{s,o) = o (103) 



where the coordinates of the events are given by 



S = (x°,x^) (104) 
O = (x^,0, a = 1,2, 3 (105) 

We compute the frequency measured by the observer, in terms of the frequency transmitted by the satellite. Prior 
to the measurement, the satellite and observer are assumed to have calibrated their equipment, when they were on a 
common world line (before launch), as discussed above. We define 

UJo — frequency measured by observer at reception event O (106) 
ujs — satellite proper frequency transmitted at event S, (107) 
as measured with respect to frequency standard on board satellite (108) 

During a proper time dr — dso/c measured by the observer, there are N cycles of the radiation received. These cycles 
were sent by the satellite during a proper time dsg/c , as measured by the satellite clock. Assuming that the signal 
does not pile up anywhere, e.g., a static space-time, we can say that 

N = ujo—=uj,'^ (109) 
c c 

Note that dso is measured at event O and that dsg is measured at event S and that these two events are arbitrarily 
separated (but connect by a null geodesic) . We can then relate the world lines of the observer and satellite by 

dsl = -g,jiO)dxldxi (110) 
dsl = -gij{S)dx\dx{ (111) 

where dx* and dx\ are increments along the observer and satellite world lines, and gij{0) and gij{S) is the metric 
evaluated at events O and S, respectively. The ratio of frequencies can then be written as 



UJo _ dss 

LUs dSo 



/ o\ dx], dxl 
9ijW) 'dX'dX 
, ^\ dx' dxi 



1/2 

(112) 



where A is a common parameter for the two world lines. Note that the frequency ratio in Eq. (|112|l depends on two 
space-time events, S and O, so this ratio is a two-point tensor, see Section V and Section X-C. 

We compute the observer's world line by taking the observer to be stationary on the surface of the rotating Earth. 
Using ECEF (Earth-centered Earth fixed) coordinates, ?/*, i = 0,1,2,3, the spatial coordinates of the observer, y°', 
are constant. The transformation from ECI Cartesian-like space-time coordinates x^, to ECEF Cartesian space-time 
coordinates y*, is of the form given in Eq. (|70|l . with uj il, where 51 is the angular frequency of rotation of the Earth, 
with respect to an ECI coordinate system with a common z-axis. In the ECI coordinates x% taking the observer to 
be located on the j/^-axis on the equator, the observer's world line is then given by choosing — R, — 0, y^ = 0, 
where R is Earth's equatorial radius. 

As a metric in the vicinity of the Earth, we choose Eq. H68|l . with x^ x° and we neglect the Earth's quadrupole 
moment, taking J2 — 0, so that V{r, 9) ^ (p — —GM/r. The world function for this metric is given by: 

n{x\,xi) = ^r,,,Ax'Ax^ + h^^Ax'Ax^'^ i^^!^) ^^^^^ 

where cAt = X2 — x^, and 0i and 62 are defined by 



a = 1,2 



(114) 
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We choose the parameter A = in Eq. H112() . so the ratio of frequencies becomes 

1/2 



1 + 5, - (1 - 5.)<5a/3f>f 



1 + So - {I - 6o)Sc,l3V^vS 



where 



[4*1 




2GM' 




o 


c^r 






2GM' 




s 


c^r 



da;" 

and the subscripts o and s on the square brackets indicate evaluation at events O and S*, respectively. 

For typical Earth orbiting satellite applications, the (S-terms and velocity terms in Eq. (|115|l are all small: 
6,r^vlr^v1r^ 0(2) whcrc 0(2) - 10-1°. 

It is then convenient to look at the fractional frequency shift, making the expansion 



(115) 



(116) 

(117) 

(118) 
(119) 

(120) 
(121) 

5o - 



— - 1 



1 



6o + vl-vl)-\{6, 



-[{5o-vlf-5oi 



(122) 



where 0(5) ^ 10^^^ for typical Earth-orbiting satellites. 

As in Eq. H78|l - (|99|l . we use the classical mechanical conditions to integrate the geodesic equations, for the current 
metric. So the satellite world line is known for the given initial conditions. The satellite world line is parametrized by 
coordinate time: x'^{x^). The observer's world line, x^{x^), is also parametrized by coordinate time. We know the 
observer's world line because we assume that the observer is located on the surface of the Earth on the axis. The 
calculation is then done as follows. For a given (emission event) satellite coordinate time, x", and spatial position, 
x", we compute the observer's reception event coordinate time x" by solving: 



f}(x°,x"(x°),xO,xf(x°)) = 



(123) 



The emission and reception event coordinates, in Eq. (|104|) and H105|) . are then known. These coordinates are used 
in Eq. (|122|l to plot the observed frequency shift, ^ — 1 as a function of the reception event time x". 

A computer program in the Mathematica programming language was written to carry out the required compu- 
tations. The program does not include the (relatively large) atmospheric signal propagation delays, although these 
effects could be included in the future work. Furthermore, no allowance is made for the Earth obscuring the signal 
from the satellite to observer. 



F. Doppler plus Gravitational Frequency Shift: Numerical Results for Various Orbits 

Proper time on board a satellite is not a linear function of coordinate time. This nonlinear functional dependence 
can be seen on a plot of the fractional frequency shift, Aw/ws, where Aw = Wo — Ws, and lOq is the satellite frequency 
as observed with respect to a reference oscillator on the geoid, and lJs is the proper frequency emitted by the satellite- 
as determined by a frequency standard on board the satellite. See Eq. (|122(l for the definition of Aw/oj in terms 
of orbital parameters. (Here we assume that the satellite frequency standard has been calibrated and that it has 
not been altered, as has been done in GPS satellites by the "factory offset".) The plots in Figures [TUlini show the 
fractional frequency shift, Aw/wc, vs. At, where At is the elapsed coordinate time measured from time of perigee, for 
the LEO, GEO, HEO, and GPS satellites. (At perigee we have arbitrarily taken At = 0.) 

The frequency shifts, similar to those shown in Figures [TriHl4l would be observed at other locations than the Earth's 
geoid. For example, whenever the proper time is nearly equal to coordinate time, cc", approximately the same plot 
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FIG. 10: The frequency shift (Doppler plus gravitational) as observed from the geoid for the LEO satellite is shown vs. 
coordinate time At = Ax^ /c in units of hours. 



would be observed. This is the case for observers in laboratories located at altitude above the geoid. Similar, but more 
complicated frequency shifts than plotted in Figures [11711141 would be observed by satellites. Note that the observed 
frequency shift depends on position and velocity of both observer and satellite. The dependence on spatial position 
of observer and satellite is due to the fact that space is not homogeneous in a curved space-time, due to the presence 
of the gravitational field. 

The LEO satellite in Figure^] has an overall shift to lower frequency, due predominantly to its high orbital speed. 
The effect of time dilation dominates the effect due to the gravitational potential difference. Superimposed on the 
overall negative frequency shift, is a small-amplitude rapid variation in the frequency shift due to the satellite's short 
orbital period. This rapid variation is due to the slight orbital eccentricity and the satellite's orbital inclination 
interacting with the Earth's gravitational field. 

The frequency shift of the GEO satellite has an overall positive value of 5.38775x10^^°, and is dominated by 
the gravitational potential frequency shift, causing a shift to higher frequency when observed from the ground, see 
Figure ITTl The plot shows the residual frequency shift, after subtraction of 5.38775x10"^", resulting from a slightly 
non-circular orbit. The satellite parameters have been chosen so that the orbital inclination is zero. The small 
variation in frequency shift is the result of choosing an initial position and velocity that do not produce a perfectly 
circular orbit (as may happen when real a satellite is inserted into orbit). In the initial data used here, the result is 
a small variation in the frequency shift on the order of 2 x 10^^*. 

In contrast, the HEO satellite in the Molniya orbit shows a markedly different behavior, see Figure^l The frequency 
shift is negative for short periods of time when the satellite is near perigee (near At = 0) but is mostly positive for 
a long time (for approximately 12 hours, roughly between At — 0.5 and At = 11.5 hours) when the satellite spends 
a long time at high altitude. The HEO satellite makes a rapid transition from the extremes of frequency shift of 
Aoj/u! — +4: X 10^^° to —4 X 10~^°, in a time on the order one hour, whereas the satellite period is approximately 
12.4 hours. 

The frequency shift for the GPS satellite is shown in Figure IT^ The GPS satellite is in a high, almost circular 
orbit, with an inclination of 55.3°. The net frequency shift is observed to be positive (higher frequencies are observed 
on the geoid) with a slight variation due to orbital eccentricity. The net shift is approximately 4.4x10""'^'^, and in 
the actual GPS system, it is compensated by the factory offset of -4.4x 10~^°. In addition to this constant frequency 
shift, the frequency shift also varies at the period of the GPS satellite. Furthermore, plots of frequency shift for all 
satellites typically have non-zero amplitudes at higher Fourier components than the satellite period. As an example, 
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FIG. 14: The frequency shifts (Doppler plus gravitational) as observed from the geoid is shown vs. coordinate time At — Ax^/c 
in units of hours, for the LEO, GEO, HEO (Molniya), and GPS satellites. 
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FIG. 15: A plot of Eq. 1125^ shows that the next higher harmonic, is present in the fractional frequency shift observed on the 
geoid for the GPS satellite. 

for the GPS case, we can fit the function 

A + Bsm{2TTft + (j)) (124) 

where A, B, lo' and (j) are constants, to be fit to the computed function Aw/ws. The values of the best fit parameters 
are given by A = 4.4736 x lO-^", B = -5.78749 x IQ-^^^ / = 1.99009 day-\ and (/)/27r = 0.25011 day. To demonstrate 
the presence of higher harmonics. Figure [T^ shows a plot of the difference 

— - (A + Bsin(27r/i + 0)) (125) 

The plot shows that the next harmonic present is at twice the basic GPS period. From Figure we see that the 
amplitude for this next harmonic, which is on the order of 1x10""'^^, is significantly smaller than for the fundamental 
frequency. Even though we have chosen a spherical gravity model for the Earth for the frequency shift calculation, 
higher harmonics are present here because the GPS orbit is slightly eccentric. In general, higher harmonics are present 
in all satellites because the orbit is eccentric and because the Earth has a non-spherical gravitational field. In the 
frequency shift calculation presented here, we have neglected the effect of the Earth quadrupole moment, however, it 
can easily be incorporated into the calculation. 



IX. GEOLOCATION IN CURVED SPACE-TIME 



Location of a source of electromagnetic radiation by using multiple receivers is termed geolocation 




The source (called the emitter) can radiate a continuous wave (cw) signal or a pulse. In either case, there are two 
methods commonly employed to locate the emitter: time difference of arrival (TDOA) and frequency difference of 
arrival (FDOA). The TDOA technique is based on differences of time of flight of the signal from the emitter to each 
receiver. FDOA is based on the Doppler effect and the difference of frequencies observed by each receiver. 

Often the emitter of interest is located on the surface of the Earth. This is a helpful mathematical constraint that 
reduces the number or receivers needed to locate the emitter. Since satellites are often used to receive signals, and 
they are costly resources, this constraint helps reduce the number of satellites required for geolocation. 
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Electromagnetic waves propagate at an (almost) constant velocity in an ECI system of coordinates. However, in 
this system of coordinates, an emitter that is stationary on the Earth's surface, is rotating with respect to the ECI 
coordinates, and hence the constraint that the emitter is on the Earth surface is actually a time-dependent constraint 
(for an Earth model that is not a surface of revolution). 



A. Time Difference of Arrival (TDOA) Geolocation 



In this section, the world function formalism is used to write the general equations for TDOA geolocation, in 
an arbitrary system of coordinates, either ECEF or ECI coordinates. Since the world function of space-time is an 
invariant, the general equations for geolocation are covariant: they are valid in any system of coordinates. I assume 
that the emitter is located at a space-time event whose coordinates are x^, i = 0, 1, 2, 3. Similarly, I assume that three 
satellites, s = 1,2,3, receive the emitter signal at space-time events with coordinates x*. The satellites' ephemerides 
are assumed known, and each satellite carries a clock so the three reception times x^, s — 1,2,3, are known as well 
as the spatial coordinates of the reception events. So all satellite space-time coordinates x]. are known. The time of 
emission at the emitter, a;^, and the spatial position of the emitter, x" , a = 1, 2, 3, constitute four unknowns. 

The electromagnetic waves from emitter to each satellite travel on null geodesies that connect the emission event 
to the three reception events, at each of the three satellites. In an arbitrary system of coordinates, the emitter 
coordinates are related to the satellite coordinates by the three equations 

n{xi,xl)=0 s = 1,2,3 (126) 

where fl is the world function of the space-time. The causality conditions requiring that the reception occur after 
emission, < x^, must also be added for a unique solution. 

The constraint that the emitter is on the Earth surface means that the emitter coordinates lie on a 4-dimensional 
hypersurface, i.e., the two dimensional surface of the Earth sweeps out a 3-dimensional hypersurface (over time) given 
by: 

xiO - (127) 

Equations (|126|l and (|127|l form a system of four nonlinear equations for the four space-time coordinates a;^ of the 
emitter. Since this system of equations is nonlinear, it may have multiple solutions, or, no solutions at all. 

If the emitter is not on the surface of the Earth, then four satellites (not three) are needed to provide four equations 
(s=l, 2,3,4) of the form of Eqs. fT^ . 

As an example of the application of TDOA to locate an emitter, assume the emitter is on the Earth surface 
and that the space-time surrounding Earth is well-modelled by the Schwarzschild metric. The world function for a 
Schwarzschild space-time is known |^| . Using Cartesian-like coordinates 1] and solving Eq. 1)1261) by iteration 
leads to 



GM 




cos 6o — cos 6 s 



^(Xs,Xo) 



(128) 



where x^ and Xo are the spatial coordinates of satellites and emitter, and |xs — Xo| is the three-dimensional Euclidean 
distance between satellite s at reception time and emitter at emission time. When three satellites, s = 1,2,3, make 
time of arrival measurements, we can form two independent TDOA equations: 

= e(xi,x„)-e(x2,x„) (129) 
= e(x2,Xo)-^(x3,Xo) (130) 

Equations H129|) - (|130|l are the TDOA equations which contain the effects of the gravitational field on delay of signal 
propagation (Shapiro time delay, which is on the order of 40 ps). To these two equations, I must add the constraint 
that the emitter is on the Earth surface. Equations H129() - (|130|l arc in Schwarschild coordinates, which are essentially 
ECI coordinates. The transformation from these ECI coordinates, x ~ {x^,x'^,x^), to rotating ECEF coordinates, 
y = {y^,y'^,y^), is given by 
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(131) 



39 



In three-dimensional notation, there exists a rotation matrix R, given by Eq. H131(l . that relates the spatial coordinates 
at a given coordinate time 

y = Bix") ■ X (132) 

The coordinate time cc" is taken to be the same in ECI and ECEF coordinates. 

For an emission event time x'^, the rotating ECEF coordinates are related to the Schwarzschild coordinates by 

Yo = R{x°) ■ xo (133) 

In rotating ECEF coordinates, the Earth's surface is given by 

/(y) - (134) 

The constraint in Eq. (|127() that the emitter is located on the Earth surface is then time-independent and can be 
written in ECEF rotating coordinates as 

f{yo) = f{R{x"J-^o)^0 (135) 

Note that the unknown time of emission x'^ enters explicitly in Eq. H135|l . This time can be eliminated by using 
Eq. 1)1281) (with s = 1) in Ea. H135|l . leading to the constraint equation 

/(i?(a;°-e(xi,x„)).x„)-0 (136) 

The three equations to be solved for the three spatial ECI coordinates Xq of the emitter are Eq. ()129)I - H13()|I and 
(|13mi . When the coordinates Xq are found, they are to be substituted into Eq. 1)1281) to compute the emission time 
Xq, in terms of the known satellite reception times x^. This time is then used in Eq. (|133|l to compute the emitter 
ECEF rotating coodinates yo- 



B. Doppler Effect in a Gravitational Field 



The previous section described the equations needed to carry out gelocation by the technique of TDOA. An 
alternative and complimentary method is based on measurements of frequency difference of arrival (FDOA). 

The frequency that an observer sees emitted by a distant source is not the same as the frequency transmitted 
by that source. More specifically, the observed frequency differs from the proper (transmitted) frequency because 
of relative motion between source and observer (Doppler effect) and because of a gravitational potential difference 
between the position of the source and position of the observer. Most often, the Doppler effect and the gravitational 
potential effect are treated as separate effects. In a static space-time, there is no difference. However, in a general 
space-time (not static or stationary), these two effects are inseparable. One expression describes both effects. Below, 
we follow closely the derivation by Synge |3ll | . Consider a space-time with an arbitrary gravitational field (not static 
or stationary) . The source of electromagnetic radiation (emitter) travels on world line C" and the observer (satellite 
with receiver) has world line C, see Figure 1161 The emitted photon has momentum in the comoving coordinate 
system of the emitter. In the comoving system of coordinates of the observer (satellite) the photon has momentum 
p*. For the comoving system of coordinates of the observer, we define a tetrad A^^^ a = 0, 1, 2, 3. The three spatial 
components of the tetrad play the role of orthonormal laboratory unit vectors. This tetrad frame does not have 
to be a Fermi frame (it does not need to be a non-rotating frame). We choose the time- like basis of the tetrad to 
be the observer 4- velocity A^g^ = V^. We also use the fact that the three space- like basis vectors of the tetrad are 

orthogonal to this timelike basis vector A^^-j Ai(o) = 0, for a — 1,2,3, and orthonormal to each other A'^^-j Xi(p) — 0, 
for a, P ~ 1,2,3 and a ^ (3. The inner product between the tetrad basis is the Minkowski matrix: 

^(a) ^t(b) = Kab) (137) 

where 77(a6) = diag(— 1, +1, +1, +1). As usual, tetrad indices are raised and lowered with rj(^ab)- 

Define the 3- velocity of the source on C, relative to observer on C, by the three invariant projections on the tetrad 
at P: 



(138) 
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source observer 

FIG. 16: The world lines of the emitter C' and receiver C are shown, together with their four velocity, at emission event 
P' , and V at reception event P. The photon at emission time has proper momentum and at reception time, has momentum 



However, the velocity of the source and the velocity of the observer cannot be compared directly because they 
are at different space-time points, P' and P, respectively. Therefore, the 4-velocity vector at P' must be parallel 
translated to P, defining the vector Vi at P by 

V, = g,,, V' (139) 

where giji is the parallel propagator from P' to P. Taking v^Vi = —1, so that is the velocity in units of c, the 
invariant components on the tetrad basis are related by 

i;(0) = (l + V^'^^^^^) ^'^ = -Vf^o) (140) 
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Then if f = 1, the other components v^"'' = 0, a = 1, 2, 3, so that, as compared at point P, the velocity of source 
and observer are equaL 

Now define at point P the unit vector orthogonal to the observer 4-velocity 

nV' = (141) 

such that the vector lies in the 2-element which contains the tangent vector to C and the tangent vector pi to 
the null geodesic P' P that connects the source and observer. Then by Eq. (|141|l . J'iA^Q^ = r(o) = 0. Define the speed 
of recession of the source with respect to the observer by 

«fl = w,r* (142) 

since r'^^ — 0. 

Since the curve P' P is a null geodesic, the emitted photon proper momentum pii is related to the photon momentum 
Pi at P by parallel transport. Furthermore, the scalar product of two 4- vectors is invariant under parallel transport: 

p,,V''^py (143) 

The energy of the emitted photon in the comoving frame of reference of the emitter is 

-E'= p'y = p(o)z;(°) + (144) 

The energy of a photon, with momentum p*, seen in the comoving frame of the observer with 4-velocity V^, is given 

by 

E = -p,V' = -p, A|o) = -p(o) = (145) 

From the definition of unit vector r% i.e., that it exists in the 2-element that contains the tangent at P to 
C, and the null geodesic tangent vector pi, and the fact that pi is a null vector, the vector pi can be written as a 
decomposition 

pi = a{V' - r') (146) 

where r^r' — 1 and the magnitude a is to be determined. Multiplying Eq. (|146|) hy Vi, I find that p'-Vi = —a = —E. 
Therefore, Eq. (|146|l can be written as 

p' = E{V'-r') (147) 

From Eq. (|144|l . and using p(o) = ^E, 

= £;w(o) (148) 

The term = A(a)^^'"' and using Eq. (|147|l for pi leads to 

= i?(F,-r,)A|„)«(") (149) 
= E [V.XI^^ - (150) 
= -£;r(„)t.(") = -EvR (151) 

since Vi^^a) ~ ^ because the frame is orthogonal, and we used Eq. (|142|) for the definition of the recessional velocity 
vr. Now use Eq. (|151|l in Eq. H148|) for P(q,)w''"' and Eq. (|140|l for w'") to obtain 

E = (152) 



where 



= ,5„^i;(")z;(« (153) 




and a, /3 = 1, 2, 3. In Eq. 1152(1 . we have the foUowing definitions: 
i?=photon energy observed by sateUite at P 

i?'=photon energy emitted in (proper) comoving frame of emitter at P' 

?;^=square of speed in terms of components on tetrad, defined by Eq. H153|l 

VR = recessional velocity defined by Eq. H142|l 

Equation (|152|l shows that the frequency of an electromagnetic signal observed by a satellite, E, is related to the 
frequency emitted, E' , by motional (velocity) effects. Of course, the role of the gravitational field in producing a 
frequency shift (gravitational red shift effect) is contained in Eq. (|152|l in the geometry of the curved space-time. 
The significance of Eq. H152|) is that it is a general expression, that is valid in an arbitrary space-time and does not 
assume a static or stationary space-time. However, Eq. (|152|) does assume that geometrical optics is valid because 
the geodesic law for light propagation was used. 



C. Observed Doppler Shifts 



The observed energy of a photon, E, is related to the observed angular frequency lo, hy E — Tilj. In a system of 
coordinates where the photon 4- momentum is pi, and the observer 4- velocity is , the observed photon frequency is 
given hy E = hoj = —piV^. Similarly, the energy of the photon emitted by the emitter, in the comoving frame of the 
emitter, is related to the frequency by E' = Tilj' = —pi'V' , where pi' is the 4- momentum at the emission event and 
y is the 4-velocity of the emitter. 

The fractional shift in frequency, between the emitted and received frequencies, due to motional doppler effects and 
gravitation is given by 

E' Pi'V^ 

The right side of Eq. H154|l can be written in terms of the derivative of world function of the space-time. The 
world function connects two points, P' and P, by a geodesic. The geodesic connecting P' and P may be expressed 
parametrically as x'(u), for z = 0, 1, 2, 3, for uq < u < ui, where P : x^{uo) and P : a:*(ui). 

The world function is a two-point scalar, depending on points P' and P. The covariant derivatives of the world 
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function with respect to point P' and P are given by 

a- = n,,{P',P) = — ^^{ui-Uo)U^, (155) 

f],; = ^,{P', P)^ — = {ui- Uo)U, (156) 

The photon momentum, pii at point P', is related to the photon momentum, pi at point P, by paraUel transport. 
Also, the photon momentum is in the direction of the tangent to the geodesic. Also, the magnitude of a vector (photon 
momentum) is constant under parallel transport. Therefore, the momentum is proportional to the tangents at point 
P' and at P 

U^' = ap.,: (157) 
U., = ap^ (158) 

where a is the same proportionality constant. Using Ea. H155| l- H158f) . the frequency difference in Eq. H154|l can be 
written as 

Eq. H159|) gives the relation between the frequency uj observed by a satellite at the reception event P and the frequency 
of the emitted electromagnetic signal at the emission event P'. The satellite has 4- velocity and the emitter has 
4-velocity . The frequency uj' is the proper frequency of the electromagnetic signal broadcast by the emitter in 
its comoving frame, i.e., the frame in which the emitter is at rest. Eq. 1)1591) contains the effect of the gravitational 
potential (gravitational red shift) as well as the effect associated with the relative motion of the source and observer, 
which is usually called the Doppler effect. 

The significance of the expression in Eq. (|159|l is that it can be evaluated in any system of coordinates (any reference 
frame). For example, Eq. (|159|l can be evaluated in inertial or rotating coordinates. The key ingredient to carrying 
out the calculation is that the world function must be computed for the given space-time. This has already been 
done for the Schwarzschild space-time, which models the Earth as a sphere [lllsJi. A similar calculation of the world 
function should be done to include the effects of the Earth's quadrupole potential J2. However, the gravitational 
effects due to J2 are three orders of magnitude smaller than the monopole terms, and for some applications, may be 
negligible when dealing with computations of signal propagation time or frequency shift. For example, if the time 
delay due to the monopole contribution is of the order of 40 ps, then the effect of J2 is expected to be three orders 
of magnitude smaller. Furthermore, the doppler effect is not a cumulative effect such as time dilation, so these terms 
do not increase in size with time. Therefore, for many purposes the world function for the Schwarzschild space-time 
is sufficient. 



D. Frequency Difference of Arrival (FDOA) Geolocation 

In this section, we derive the basic relations for geolocation by frequency difference of arrival (FDOA) in a curved 
space-time. Geolocation means that measurement of the frequency of an emitter by several satellites (in relative 
motion, and at known positions in space-time) can be used to locate the emitter in space-time. Here we will take into 
account the frequency shifts due to the relative motion of satellite and emitter, as well as due to the gravitational 
potential effect. When satellites in different orbital regimes, e.g., LEO, GEO and HEO, are combined, the gravitational 
potential differences can lead to significant errors in emitter positions. We neglect the (relatively large) effects that 
can result from atmospheric propagation delays. 

As derived in section X, subsection G, Eq. (|159|l gives the basic expression that connects the frequency emitted 
(proper frequency) with the frequency observed, taking account the relative motion as well as different gravitational 
potentials. Assume that the emitter sends a signal at event Pq — {x^, x") and that this signal is received by a satellite 
at event P, = (a;°,a;"). The emitted frequency as measured at the emitter using a reference oscillator is lUq- The 
frequency observed at the satellite, lUs, is observed by using a reference oscillator at the satellite identical to that used 
at the emitter to determine the emitter's frequency. The fractional frequency difference can be written in terms of 
the world function and the 4-velocities of the emitter and satellite receiver: 
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where for practical applications i? << 1. In Eq. (|160|) . we have the following definitions 



ft,. 



yi„ _ _ 4_velocity of emitter 
ds 

dx^ 

yi, ^ ^ 4-velocity of satellite 
ds 

ujo — proper frequency of emitter, with 4-veIocity V^*° 
LUs — frequency of emitter, as observed at 
satellite with 4-velocity V^" 

The frequency observed by a given satellite can be written as 

iv,=Uo{l-R) (161) 

where R depends on the 4- velocities of emitter and receiver, and the world function of space-time. For the Minkowski 
metric given in Appendix B, explicit evaluation of Eq. (|16U|I gives 

7s(l-n-v,) 

= Wo = — (162 

7o(l - n - Vo) 

where Vq and Vg are the velocities of the emitter and satellite (in units of v /c) , and the unit n vector is defined as 
n = [ys — Vo) /\ys — Vo\, where Tq and are the spatial position of emitter and satellite at points of emission and 
reception. 

The measured quantities are frequency differences, which can be formed from difference of satellite frequency given 
in Eqs. H161|l . We do not give the final complicated expressions here for frequency differences. However, oscillator 
frequency errors clearly contribute to errors in estimated emitter position. For example, preliminary analysis shows 
that the order of magnitude error in emitter location is given by 

5uj |x^-Xo| 

(5xo = — 7j- (163) 

where Xq, Vo and x^, are the position and speed of the emitter and receiver, respectively, and ^ is the fractional 
frequency error. For example, for the typical values of satellite speed Vs ~ 10~^, and satellite oscillator frequency 
error ^ ~ 10"^*^, the order of magnitude of emitter position error is |(5xo| ~100 m. A frequency error on the order 
of 10"^" can also be expected to occur when a LEO and a GEO satellite are used together for geolocation without 
compensating for the gravitational frequency shift. In such as case, the use of two satellites can lead to a similar 
position error on the order of 100 m. Detailed investigation of Eq. H161|l for high-accuracy geolocation is left for future 
work. 



X. SUMMARY 



In this article we have considered the elements of the general problem of navigation in space-time, as well as the 
restricted problem of clock synchronization, within the context of a metric theory of gravity, such as general relativity. 
In most real applications, such as the GPS, a user is interested in determining his space and time coordinates, rather 
than just time. General relativity deals with the effect of motion and gravitational potential differences on clocks, and 
it highlights features of the space-time navigation problem that must be included in any (classical or quantum) theory 
of navigation (or clock s y nchronization). Several quantum mechanical schemes have been proposed to synchronize 
clocks m in m Hi, mm |2g |33. At the present time, the effects of motion and gravitational potential 
differences have not been explicitly incorporated into these quantum approaches to clock synchronization. Clearly the 
magnitude of the relativistic effects is such that it must be considered in future quantum approaches to navigation. 
At the present, there is no quantum theory of navigation in space-time (analogous to the classical theory in Ref. j^) 
which permits a user to determine their four space-time coordinates. 
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FIG. 18: The time and position of an electromagnetic emitter is shown at a space-time point {t,r), together with satellites, 
Sa, located at {Ta, Ra), a = 1, 2, 3, 4. 
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APPENDIX A: SIGNIFICANCE OF ACCURATE CLOCK SYNCHRONIZATION FOR GEOLOCATION 

Consider the problem of locating an emitter of electromagnetic radiation located at a space-time point (t, r) near 
the Earth's surface. The time of emission and position of this emitter can be determined from signals received at 
fours satelhtes, located at reception events (Ta,Ra), a — 1,2,3,4, see Fig.lTHl 

In an inertial system of coordinates (such as the ECI frame), the time and position of the emitter can be computed 
by solving the four equations 

c(r,-t) = |r-R„|, a =1,2,3,4 (Al) 

for {t, r), where the satellite positions, Ra, and reception times, Ta, are known, and c is the speed of light in vacuum. 
In this Appendix, for simplicity, we neglect the distinction between proper time and coordinate time. 

By taking cyclic differences of Eqns. (|A1(I . we cancel out the emission event time t and obtain time difference of 
arrival (TDOA) equations 

c{Ti-T2) = |r-Ri|-|r-R2| (A2) 
c{T2'T^) = |r-R2|-|r-R3| (A3) 
c{T2-Ti) = |r-R3|-|r-R4| (A4) 

where the measured quantities are the time differences Ta — Ta+i at the satellites. The Eqns. I|A2|I can be solved for 
the coordinates of the emitter r, in an inertial coordinate system. The value of r must then be substituted into one 
Eqns. (|A1|I to obtain the emission time t. The coordinates of the emitter, r, can then be transformed to geocentric 
(rotating with the Earth) coordinates, r', by a time-dependent rotation 

r' = D(t) • r (A5) 

where D(i) is the transformation matrix. 

A similar problem occurs if we know that the emitter is located on the surface of the Earth, so that we have the 
constraining equation 



|r| = TE 



(A6) 




FIG. 19: The orientation of an electromagnetic emitter O is shown with respect to two satellites 5*1 and 5*2. 



where r£;=6378 km is the Earth's radius. In this case, a minimum of three satelhtes are sufficient to locate the emitter. 
The first two of Eqns. HA2|I are solved together with the constraining Eqns. (|A6() for r. Then the geocentric (rotating 
with the Earth) coordinates, r', are found by the time-dependent rotation in Eqns. (|A5p . 

The question then arises: How does the accuracy of the synchronization of the clocks affect the accuracy of 
geolocation? Assume that the clocks of two satellites are synchronized, and that the satellite positions are known. 
The maximum TDOA occurs for satellites separated by an angle (j) given by (see Figure IHIjl 

sincf) = 12/rs (A7) 

The path length difference between these two satellites is 



M = l2-li = ^rl - r| - (r, - r^) (A8) 
The maximum time difference of arrival (TDOA) for a signal is 



^^^A/^^^^ Vrg-4-(r.-r^ (A9) 
c c c 

For geosynchronous satellites, where =6378 km and rs =42,164 km, the maximum (TDOA) is At =19.64 ms, and 
it occurs when when one satellite has a longitude that is the same as the emitter and the other satellite is at (j) =81.3°, 
see Figure El Due to geometric constraints, the whole range difference is contained in the maximum time delay of 
19.64 ms. If the clocks in the two satellites are synchronized only to an accuracy of, say 10 ns, then the whole range 
difference is compressed into a time delay of 19.64 ms. The fractional range error, and consequently the order of 
magnitude of the position error (neglecting geometric dilution of precision (GDOP) factors) is given by 

Ax - A/ = (I2 - h)-^^ 3000m = 1.5nm (AlO) 

19.64ms 



An improvement in the clock synchronization, say by three orders of magnitude, translates directly into three orders 
of magnitude improvement in position accuracy, resulting in position accuracy of the order of 3 m. 
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In the discussion in this appendix, we have neglected atmospheric (tropospheric and ionospheric) time-delay effects, 
accuracy of computational correlation algorithms, and noise in the receiving system. An implicit assumption is also 
made that the system of coordinates (inertial reference frame) is accurate, and that the time-dependent transformation 
between inertial and geocentric (rotating) coordinates is accurately known. 

APPENDIX B: CONVENTIONS AND NOTATION 

Where not explicitly stated otherwise, we use the convention that Roman indices, such as found on space-time 
coordinates x% take the values i = 0, 1, 2, 3, and Greek indices take values a = 1, 2, 3. Summation is implied over the 
range of any index when the same index appears in a lower and upper position. In some cases, summation over Greek 
indices is implied when indices both appear in upper position, such as in x°'dx°'. 

If and -I- dx* are two events along the world line of an ideal clock, then the square of the proper time 
interval between these events is dr — ds/c, where the measure ds is given in terms of the space-time metric as 
ds^ = —gij dx^ dxK I choose to have the signature -1-2. When gij is diagonalized at any given space-time point, 
the elements can take the form of the Minkowski metric given by Tyoo = "li Vap = Sap- In discussion of an observer 
carrying a tetrad, the Mnkowski matrix, ri(^ab) — Vab, is used. 
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